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ABSTRACT 
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Edge  extraction  is  an  image  processing  technique  for  defining 
the  edge  information  in  an  image.  This  effort  researches  different 
edging  processes  as  applied  to  preprocessing  for  two  pattern 
recognition  processes.  The  first  one  is  a  cross-correlation  method 
to  find  a  target  given  that  the  target' has  a  known  size,  orientation, 
and  aspect.  Correlation  is  performed  in  the  spatial  frequency  domain 
with  two-dimensional  fast  Fourier  transforms  of  the  searched  edge 
image  and  a  hand  drawn  edge  template  to  correct  for  translation  only. 

The  second  pattern  recognition  process  researched  also  uses 
edging  as  one  step  of  a  purely  spatial  domain  algorithm.  The 

approach  locates  targets  in  infrared  images  that  can  be  described 

-J? 

a ^hot^  clusters.  A  cluster  recognition  algorithm  by  Hamadani  is 
implemented  and  altered  for  testing  of  local  thresholding  and 
thresholding  rules.  The  algorithm  is  shown  to  be  effective  on 
real  infrared  images  provided  by  the  thesis  sponsor,  the  Air  Force 
Armament  Laboratory  %,  Eg  1  in  AFB,  Florida. 


IMAGE  FILTERING  WITH 
BOOLEAN  AND  STATISTICAL  OPERATORS 

I.  INTRODUCTION 

BACKGROUND 

An  image  is  the  representation  of  the  electro-magnetic  energy  at 
a  certain  time  in  a  plane  in  space.  Typically  the  energy  is  in  the 
spectrum  of  visible  light  or  infrared  light.  Image  pattern  recogni¬ 
tion  is  the  process  of  analyzing  an  image  to  find  known  patterns 
that  aid  in  the  classification  of  the  image  information  content. 

Image  pattern  recognition  has  many  applications  such  as:  (Ref  1) 

1)  document  processing 

2)  industrial  automation 

3)  medicine  and  biology 

4)  military  target  finding,  guidance. 

There  is  no  pattern  recognition  machine  to  date  that  approaches 
human  capabilities.  Pattern  recognition  problems  are  difficult 
because  of  the  natural  complexity  of  the  world.  Humans  use  form 
(gestalt)  matching  whereas  all  machines  perform  template  matching 
when  the  form  has  a  specific  rendition.  Machines  can  work  practically 
under  controlled  conditions.  For  example,  retail  stores  often  use 
bar  code  labels  on  their  products  for  inventory  control.  Lasers 
can  read  these  carefully  printed  labels  if  they  are  held  at  the 
proper  distance. 

The  general  approach  for  recognition  of  objects  is  to  extract 
a  set  of  features  that  describe  the  object  and  search  the  scene 
for  a  match.  The  set  of  features  is  the  template  and  the  searching 


process  is  called  template  matching.  The  selection  of  features  will 
determine  the  success  of  the  match.  (Ref  2) 

The  mechanics  of  the  matching  process  must  employ  some  decision 
theory  that  will  allow  for  the  unavoidable  errors.  The  acceptance 
region  for  a  particular  feature  class  should  be  as  large  as  possible 
without  overlapping  the  region  for  a  different  class.  Often  a 
confidence  rating  will  be  assigned  to  each  match. 

Segmentation  is  the  process  of  dividing  the  image  into  meaningful 
parts  during  the  analysis.  Meaningful  parts  might  be  sections  of 
the  scene  where  there  is  a  likelihood  of  finding  the  sought  objects; 
or  they  might  be  regions  that  correspond  to  parts  of  the  sought  object 
When  a  process  has  the  ability  to  determine  how  to  beneficially 
change  its  version  of  the  template,  then  we  say  that  it  has  learning 
ability.  Learning  should  improve  the  feature  set  or  it  will  fail. 
Artificial  intelligence,  a  field  similar  to  pattern  recognition, 
lends  much  to  learning  theory. 
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SCOPE 

This  thesis  works  with  digitized  images  of  visible  light  and 
infrared  light  quantized  into  16  grey  levels.  The  video,  or  visible 
light,  images  are  represented  in  256  X  256  matrices,  while  the  infrared 
images  are  represented  in  96  X  100  matrices.  The  video  images  were 
digitized  in-house  from  photographs.  The  Air  Force  Armament 
Laboratory,  Eg! in  AF8,  Florida,  provided  the  infrared  images.  These 
images  portray  real  scenes  with  military  application. 

Images  may  be  described  by  many  features,  but  common  ones  are 
edges,  texture,  regions,  shapes,  and  spatial  frequency  content.  The 
first  four  listed  are  defined  in  the  spatial  domain,  while  the  last 
is  defined  by  a  transformation  such  as  the  Fourier  transform. 

Preprocessing  is  the  filtering  or  reduction  of  an  image  before 
the  template  matching  process.  This  thesis  is  concerned  mainly 
with  the  use  of  linear  and  nonlinear  spatial  filters  for  the  pre¬ 
processing  of  an  image.  Especially,  the  filters  will  be  designed  to 
extract  the  edges  in  the  image.  This  is  similar  to  high  frequency 
enhancement  in  the  frequency  domain. 

Some  template  matching  processes  are  considered,  and  the 
edge  extraction  is  optimized  for  each  template  matching  process. 

While  edge  extraction  is  emphasized,  no  fitting  of  mathematical 
models  for  texture  or  shapes  is  used. 

The  techniques  to  be  used  will  be  chosen  to  be  as  immune  as 
possible  to  clutter.  Clutter  is  all  of  the  background  scene  that 
distracts  the  recognition  process  from  finding  the  object.  Clutter 
may  also  be  another  object  that  occludes  the  sought  object  because 
of  their  relative  positions. 


1-3 


Present  applications  of  image  pattern  recognition  work  under 
constraints  that  depend  on  a  prior  knowledge  of  the  situation  at 
hand.  The  general  problem  does  not  assume  these  constraints.  This 
thesis  does  not  attempt  to  solve  the  complete  general  problem,  rather 
it  researches  techniques  that  are  hoped  to  contribute  to  the  solution 
of  the  general  problem. 

While  new  techniques  are  tried  on  real  scene  images,  the  results 
are  visually  inspected  for  information  content.  The  information 
content  in  question  is  in  the  edges  of  the  image.  Visual  inspection 
and  careful  observation  replace  any  mathematical  development  to 
explain  what  exactly  the  process  did.  Hypotheses  of  beneficial 
operations  are  developed  without  regard  to  linearity  or  causality. 

This  effort  assumes  the  perspective  of  the  target  is  fixed  and 
known.  Therefore,  geometrical  transformations  or  stereomapping  are 
not  considered.  Segmentation  techniques  are  also  neglected. 
Furthermore,  the  analysis  is  static,  i.e.:  successive  images  are 
not  compared.  Hence,  no  learning  ability,  movement  detection,  or 
three  dimensional  analysis  is  discussed. 


APPROACH 


Several  edge  extraction  or  edging  algorithms  are  evaluated  for 
their  representation  of  important  image  information  content.  The 
definition  of  important  information  is  subject  to  the  data  set,  and 
to  the  intended  processing  subsequent  to  the  edging  process.  In 
general,  edges  are  extracted  from  areas  with  mono ton ically  changing 
brightness  over  a  small  area  (typically  3X3  mask  neighborhood). 

The  edges  should  also  be  optimized  for  the  template  matching  algorithm 
used  in  the  recognition  process.  Modifications  to  the  evaluated 
edging  algorithms  are  tried.  The  resulting  images  are  examined 
again  for  important  image  information  content. 

Clutter  or  noise  reduction  is  in  the  form  of  thresholding.  The 
thresholding  is  typically  some  boolean  function  based  on  statistical 
operators.  The  statistics  may  be  global  or  local.  Global  thresholds 
depend  on  the  statistics  of  the  entire  image,  whereas  local 
thresholds  depend  on  a  local  neighborhood  of  the  pixel  being 
processed. 

Edged  image  values  range  from  zero  (representing  no  edge)  to 
fifteen  (representing  the  highest  contrast  edge).  When  edge  image 
pixels  are  thresholded,  the  resulting  value  is  either  zero  or  non¬ 
zero.  A  value  of  zero  precludes  further  processing  for  that  pixel. 

A  non-zero  valued  pixel  after  thresholding  can  take  on  its  original 
edge  value  or  some  function  thereof. 

The  Hamadani  algorithm  is  implemented  and  evaluated  for  techniques 
that  can  be  applied  to  more  general  pattern  recognition  problems. 

(Ref  3)  The  Hamadani  algorithm  works  solely  in  the  spatial  domain 


1-5 


for  the  purpose  of  detecting  "hot"  clusters  in  an  infrared  image. 
Modifications  to  the  Hamadani  algorithm  are  tested  to  see  if  target 
detection  can  be  improved.  Among  the  proposed  modifications  are 
improved  edging  algorithms,  local  thresholding,  and  faster  implemen¬ 
tation  techniques. 

A  more  complicated  template  matching  process  than  cluster 
detection  is  used  to  operate  on  edged  images  to  test  the  effectiveness 
of  the  edging  algorithms.  (Ref  4)  Cross-correlation  in  the 
frequency  domain  with  two-dimensional  fast  Fourier  transforms 
(2DFFT)  should  correct  template  translation.  Cromer  concludes  that 
scene  brightness  cannot  be  used  successfully  in  this  template 
matching  process,  and  recommends  preprocessing  with  edging  routines. 

New  template  matching  processes  are  contemplated  to  more 
effectively  use  edged  image  information.  Proposed  techniques  are 
discussed  as  possible  subjects  of  more  research. 


EQUIPMENT 

All  thesis  work  is  done  in  the  Signal  Processing  Laboratory 
at  the  Air  Force  Institute  of  Technology.  All  software  is  written 
in  Fortran  5  to  be  run  on  the  Eclipse  S/250  computer  system.  The 
system  has  complete  video  input  and  output  facilities.  (Ref  5) 
Output  can  be  achieved  on  a  video  monitor  or  a  Printronix  300  line 
printer.  The  linepr inter  output  is  used  to  represent  images 
for  inclusion  in  this  thesis.  Input  from  a  vidicon  camera  and 
output  to  the  video  monitor  are  executed  on  the  Nova/Cromenco 
system  which  shares  disk  memory  with  the  Eclipse  S/250  system. 


Support  software  routines  are  included  in  Appendix  A. 


II.  EDGING  ALGORITHMS 


BOOL PASS 

At  the  start  of  this  thesis  effort,  it  was  believed  that  there 
was  some  merit  to  a  proposed  edging  routine  devised  by  Blaine  Feltmate 
(Ref  6).  The  operation,  entitled  Boolpass,  consists  of  three  basic 
steps : 

1)  average  or  low  pass  filter  the  original  image 

2)  negate  the  averaged  image 

3)  threshold  using  the  criteria  that  the  original  image  and  the 
negated  averaged  image  must  be  within  a  small  range  of  each  other 
on  the  gray  scale. 

Boolpass  is  almost  equivalent  to  thresholding  the  original 
image,  around  the  middle  of  the  gray  scale.  The  operation  does  a 
selective  thresholding  of  the  middle  of  the  gray  scale,  though. 

That  is,  only  areas  where  the  average  image  is  in  the  middle  of  the 
gray  scale  remain  after  thresholding. 

For  certain  pictures,  the  operation  tends  to  reduce  the  image 
to  important  edges.  Important  edges  are  those  edges  that  are 
critical  to  the  definition  of  the  target.  The  operation  left  many 
areas  containing  no  edges  unreduced,  however.  It  also  eliminated 
some  edges  that  were  present  in  some  areas  that  had  no  values  in 
the  middle  of  the  gray  scale. 

Many  variations  were  tried  on  Boolpass  as  suggested  by  Feltmate, 
and  some  others,  too.  Nothing  seemed  to  be  able  to  leave  all 
important  edges  and  reduce  unimportant  areas.  Then,  using  the 
criterion  that  deviation  for  the  neighborhood  around  the  pixel 


being  operated  on  must  be  over  a  certain  value,  Bool  pass  gained  the 
ability  to  reduce  areas  with  middle  scale  values  and  with  no  edges. 
However,  edges  not  in  the  middle  of  the  gray  scale  still  were  not 
present.  I  then  concluded  that  Bool  pass  was  no  longer  viable  and 
pursued  the  idea  that  deviation  was  the  feature  that  defined  edges. 

DEVIATION 

Deviation  works  well  at  defining  the  edges  (see  Figures  1  and 
2).  Specifically,  mean  deviation  was  used  instead  of  standard 
deviation  because  of  faster  computation  and  noise  immunity.  Large 
masks,  or  neighborhoods,  tend  to  blur  the  edges  (see  Figure  3). 
Therefore,  most  edge  processing  uses  3X3  masks. 

Appendix  B  contains  the  listings  of  the  Fortran  programs  used 
for  edging.  The  program  used  to  calculate  the  averages  and  mean 
deviation  on  local  neighborhoods  is  STATISTICS. 

While  deviation  works  well  at  defining  the  edges,  random 
noise  at  some  parts  of  the  image  that  contain  no  edges  is  also  extracted 
to  be  part  of  the  resulting  image.  We  do  not  want  noise  to  be 
part  of  an  edged  image. 


MASK  EDGING 


Kirsch's  operator  was  tried  because  of  recommendation  by 
Hamadani  (Ref  3).  Coefficients  were  changed  to  effectively  use  the 
gray  scale  range,  but  relative  ratios  remain  the  same  (5  to  3). 

There  are  eight  mask  orientations  to  be  convolved  with  the  neighbor¬ 
hood  of  each  pixel.  The  maximum  magnitude  of  the  eight  convolutions 
is  assigned  to  the  pixel  in  the  edged  image.  That  maximum  occurs 
for  the  orientation  that  most  closely  aligns  with  the  edge  in  the 
original  image.  The  result  of  performing  Kirsch's  operator  on 
Figure  1  is  shown  in  Figure  4.  The  Fortran  program  for  the  Kirsch 
operator  is  KEDGE,  and  the  listing  is  in  Appendix  B. 

Mask  edging  is  the  process  of  convolving  several  masks  with  the 
neighborhood  of  the  pixel  being  processed  (Ref  7:97).  The  masks 
represent  model  edge  neighborhoods  at  various  orientations.  The 
maximum  convolution  of  the  masks  becomes  the  edge  value  of  the 
output  pixel.  Notice  that  the  orientation  of  the  edge  will  be  known 
after  the  maximum  is  found. 

A  mask  edging  operator,  henceforth  referred  to  as  the  wedge 
operator,  was  defined  on  a  3  X  3  mask  with  coefficients  adding  to 
zero  so  that  the  output  will  be  zero  for  a  pixel  that  is  not  near 
an  edge.  Figure  5  shows  the  masks  for  the  wedge  operator  at  the  four 
orientations.  The  actual  coefficients  are  less  but  the  relative 
ratios  remain  the  same  (2  to  1). 

Let  us  look  at  the  masks  for  orientation  #1  of  the  wedge 
operator.  These  masks  are  models  of  horizontal  edges  where  the 
center  pixel  is: 


1)  negative  (relatively)  with  negative  values  below 

2)  negative  with  negative  values  above 

3)  positive  with  positive  values  below 

4)  positive  with  positive  values  above 

Orientations  #2  and  #4  are  opposing  diagonals  and  #3  is  vertical. 
Notice  that  half  of  the  masks  are  negatives  of  the  others. 

In  order  to  save  the  information  about  which  mask  orientation 
produced  the  maximum  convolution  value,  the  wedge  operator  separates 
the  output  into  four  different  images.  Each  output  image  is  zero 
filled  initially.  If  the  mask  with  maximum  convolution  value  is  at 
orientation  #1,  the  edge  pixel  value  will  be  assigned  to  edge  image 
#1 .  Edge  pixels  at  orientation  #2  are  assigned  to  edge  image  #2. 
Similarly,  edged  images  #3  and  #4  are  formed.  If  two  or  more  masks 
at  different  orientations  produce  the  same  maximum  convolution 
value,  the  edge  output  goes  to  each  of  their  corresponding  edge 
images. 

Other  methods  to  save  the  orientation  information  can  be 
devised  to  be  more  efficient.  The  wedge  operator  separates  the  edge 
images  for  use  in  a  scheme  for  correlation  described  later.  The 
program  for  the  wedge  operator,  WEDGE,  is  listed  in  Appendix  B.  A 
program  used  to  recombine  the  edge  images  is  NCOMB  in  Appendix  A. 
Figures  6-11  show  results  of  the  wedge  operator. 


III.  TEMPLATE  MATCHING  WITH  EDGES 


MATCHING 


To  prove  the  utility  of  edged  images  in  template  matching, 
real  scenes  were  searched  using  an  existing  cross-correlation 
scheme  (Ref  8:480).  The  template  is  in  the  form  of  an  edge  image 
of  the  object  of  interest  which  is  to  be  found  in  the  searched 
edged  image.  The  template  version  of  the  object  must  have  the  same 
size,  perspective,  and  orientation  of  the  object  embedded  in  the 
searched  scene.  The  cross-correlation  scheme  will  only  correct  for 
translation  along  the  cartesian  coordinates. 

The  cross- correlation  scheme  operates  in  the  spatial  frequency 
domain.  The  two-dimensional  fast  Fourier  transforms  (2DFFT)  of  the 
searched  image  and  the  template  image  are  formed.  The  2DFFT' s 
represent  the  spatial  frequency  content  of  the  images.  Then  the 
2DFFT's  are  point  by  point  multiplied  to  form  the  2DFFT  of  the  cor¬ 
relation  image,  which  is  inverse  transformed  to  the  spatial  version 
of  the  correlation  image.  The  programs  to  form  the  2DFFT,  inverse 
the  2DFFT,  and  perform  complex  conjugate  multiplication  are  DIRECT, 
INVERSE,  and  CMULT.  Descriptions  of  DIRECT  and  INVERSE,  and  the 
listing  of  CMULT  are  in  Appendix  A. 

Since  DIRECT,  INVERSE,  and  CMULT  require  the  images  in  complex 
format,  programs  are  needed  to  convert  between  packed  video  format 
(4  bit  integer)  and  complex  format.  VDTOCP  converts  from  packed 
video  to  complex.  CPTOVD  converts  from  complex  to  packed  video. 

In  addition,  CPTOVD  will  also  locate  the  maximum  value  of  the  image 
to  be  converted.  The  maximum  will  be  of  interest  when  the  image  is 


a  correlation  image.  The  listings  of  VDTOCP  and  CPTOVD  are  in  Appendix 
A. 

The  maximum  value  in  the  correlation  image  indicates  the  row  and 
column  translation  of  the  template  image  required  to  register  with  the 
most  closely  matching  part  of  the  searched  image.  Since  the  2DFFT 
assumes  periodicity  of  the  image,  translation  causes  wrap-around  of 
end  points.  That  is,  points  moved  down  past  the  bottom  boundary  go 
to  the  top  of  the  image,  and  points  moved  past  the  right  boundary 
go  to  the  left  boundary  of  the  image. 

More  important  than  the  maximum  value,  the  peaks  in  the  cor¬ 
relation  image  identify  areas  that  are  close  to  the  template  in 
match.  The  shape  of  the  peaks  is  important,  too.  Sharp  peaks 
indicate  a  sure  match.  Dull  peaks  indicate  a  close  but  unsure  match. 


TEMPLATE  GENERATION 

The  template  image  is  composed  of  lines  at  the  edges  of  the 
object  of  interest.  The  template  is  zero  everywhere  else.  The  lines 
are  assigned  a  maximum  value  of  fifteen.  A  variable  weighting 
method  can  be  used  to  represent  the  importance  of  an  edge  in  the 
matching  process,  but  is  beyond  the  scope  of  this  thesis.  The 
generation  of  the  template  image  is  not  a  trivial  task. 

In  order  to  generate  a  template  with  the  same  scale,  orientation, 
and  perspective  as  the  object  in  the  searched  image,  the  template 
must  be  digitized  at  the  same  time  the  searched  picture  is  digitized. 
The  templates  were  digitized  from  a  transparency  film  with  a  clean 
white  sheet  of  paper  behind  it.  The  procedure  for  digitizing  the 
searched  scene  and  the  template  follows: 

1)  Tape  the  transparency  to  the  picture  near  the  bottom  so  that 
it  can  be  folded  down  out  of  the  way  when  digitizing  the  picture. 

2)  Trace  with  a  fine  tipped  transparency  marker  the  edges  of 
the  object  a f  interest. 

3)  Set  up  the  picture  with  transparency  taped  to  it  in  front 
of  the  digitizing  camera  and  adjust  camera. 

4)  Digitize  the  picture  with  the  transparency  folded  down  out 
of  the  way. 

5)  Without  moving  anything  except  for  flipping  the  transparency 
back  up  in  front  of  the  picture,  (the  trace  should  still  be  aligned 
with  the  picture)  digitize  the  transparency  with  a  clean  sheet 

of  white  paper  behind  it. 

The  digitized  template  has  to  be  cleaned  up  because  the  lines 
will  not  have  values  of  fifteen  and  the  background  will  not  be  zero. 


The  program  TONER,  listed  in  Appendix  A,  was  used  to  map  pixel  values 
less  than  eight  to  fifteen,  and  to  map  pixel  values  greater  than  or 
equal  to  eight  to  zero.  This  mapping  negates  the  digitized  version. 
Negation  is  necessary  because  the  digitizer  assigns  a  value  of  fifteen 
for  full  brightness,  and  the  lines  on  the  transparency  are  dark. 

The  template  image  after  clean-up  should  be  a  line  trace  version 
of  the  object  with  the  same  orientation,  size,  and  location  of  the 
object.  The  template  can  then  be  moved  with  a  known  translation.  The 
program  NMOVE,  listed  in  Appendix  A,  was  used  to  translate  the 
template  to  the  center  of  the  image.  Figure  12  shows  the  finished 
template  created  for  the  image  shown  in  Figure  6. 


RESULTS 


The  image  shown  in  Figure  6  was  searched  using  the  template 
shown  in  Figure  12.  The  cross-correlation  image  is  shown  in  Figure 
13.  The  peaks  are  fairly  sharp  and  the  maximum  corresponds  to  a 
template  translation  that  is  one  pixel  location  from  the  correct 
translation. 

Figure  6  is  a  good  test  for  the  discrimination  ability  of  the 
matching  proces  since  there  are  several  other  objects  of  interest 
at  approximately  the  same  size,  orientation,  and  perspective.  The 
other  peaks  in  Figure  13  that  do  not  contain  the  maximum  correspond 
to  the  other  objects  of  interest  in  the  searched  image.  The  other 
peaks  have  maximum  values  that  are  90%  or  better  of  the  global 
maximum. 

Several  other  pictures  were  searched  similarly  and  the  targets 
were  located  within  one  or  two  pixel  locations.  These  results  are 
in  Appendix  C. 


SEPARATED  ORIENTATION  CORRELATION 


With  hopes  of  improving  the  correlation  scheme  with  edges,  the 
separated  edge  orientation  method  was  combined  with  the  cross¬ 
correlation  method.  The  method  of  separating  edge  output  into  four 
different  images  with  the  wedge  operator  is  described  by  the  section 
on  mask  edging  in  Chapter  II.  We  want  to  correlate  the  edge  image 
at  orientation  #1  of  the  searched  image  with  the  edges  at  orientation 
#1  of  the  template.  Similarly,  we  want  to  correlate  the  edges  at 
orientations  #2,  #3,  and  #4.  Then,  add  the  correlation  output  for 
each  of  the  orientations  to  form  the  final  correlation  image. 

To  separate  the  template  image  into  edges  at  four  different 
orientations,  a  line  detection  and  separation  operator  was  devised. 

The  tedge  operator,  as  hereby  named,  is  very  similar  to  the  wedge 
operator  in  that  it  uses  convolution  masks.  Each  of  the  four  masks 
shown  in  Figure  14  is  convolved  with  the  neighborhood  of  the  pixel 
being  processed.  The  mask  orientation  that  produces  the  maximum 
convolution  value  determines  which  output  image  the  pixel  will  be 
assigned  to.  If  two  or  more  masks  produce  the  same  maximum,  the 
edge  pixel  will  be  assigned  to  each  corresponding  output  image. 

The  program  for  the  tedge  operator  is  TEDGE  and  is  listed  in  Appendix 
B.  Figures  15-18  show  the  separated  edge  orientations  for  Figure  12. 

The  correlation  is  performed  as  described  in  the  previous 
section  on  matching.  The  edges  are  correlated  at  the  separate 
orientations  and  then  the  correlation  results  are  summed.  The 
program  to  sum  the  correlation  results  when  they  are  in  complex  format 
is  ADDCMP  as  listed  in  Appendix  A.  Figure  19  is  the  result  of 
performing  the  separated  edge  orientation  correlation  of  Figures 
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Figure  14  -  Convolution  Masks  for  Tedge  Operator 
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IV.  CLUSTER  RECOGNITION 


HAMADANI  ALGORITHM 

In  1981,  N.A.  Hamadani  developed  an  automatic  target  cuer  for 
infrared  images.  The  algorithm  is  basically  a  cluster  recognition 
process  that  operates  solely  in  the  spatial  domain.  The  clusters  in 
the  infrared  images  are  "hot"  blobs  that  typically  represent  man¬ 
made  objects  with  much  thermal  activity,  such  as  trucks  or  tanks. 
Hamadani 's  algorithm  was  implemented  here  in  the  Signal  Processing 
Laboratory  with  the  intention  of  researching  improved  preprocessing 
techniques. 

The  algorithm  has  seven  basic  steps  as  described  below.  The 
first  three  steps  are  considered  preprocessing.  Only  the  first 
four  steps  were  actually  implemented  with  the  assumption  that,  given 
the  understanding  of  the  remaining  steps,  the  results  of  the  entire 
algorithm  may  be  known.  The  steps  are: 

1)  Enhance  the  original  image  with  a  convolutional  mask 
operator  that  has  the  same  pattern  as  the  Kirsch  operator,  but  the 
coefficients  are  10  and  -1  Instead  of  5  and  -3  respectively. 

2)  Edge  extract  the  enhanced  image  with  the  Kirsch  operator. 

3)  Threshold  the  enhanced  image  using  a  global  conjunctive 
thresholding  method.  Global  statistics  form  the  thresholds  for 
the  enhanced  and  edged  Images.  Where  both  the  enhanced  image  is 
above  its  threshold  and  the  edged  Image  is  above  its  threshold, 
a  value  of  fifteen  Is  assigned  to  the  output.  Otherwise,  the 

output  at  that  pixel  location  is  assigned  a  value  of  zero.  Conjunc¬ 
tion  refers  to  the  thresholding  on  both  the  enhanced  and  edged  images. 


The  global  threshold  for  each  is  its  global  mean  plus  its  global 
standard  deviation. 

4)  Reduce  the  thresholded  image  using  a  connectedness  test  that 
states  a  non-zero  pixel  will  be  removed  if  it  does  not  have  at  least 
two  non- zero  neighbors  above,  below,  to  the  left,  or  to  the  right. 
Diagonal  neighbors  are  ignored.  The  test  is  repeated  until  no 
pixels  are  removed. 

5)  Segment  the  non-zero  pixels  into  target  clusters  by  finding 
which  ones  are  connected  and  assigning  them  a  target  identification 
number  to  keep  track  of  which  pixels  are  in  a  particular  cluster. 

6)  Erode  the  outside  boundaries  of  the  target  clusters  with 
what  Hamadani  calls  thinning.  For  thinning,  the  desired  target 
shape  is  assumed  to  be  rectangular.  The  top  row  of  each  cluster  is 
removed  if  the  one  below  It  has  more  non-zero  pixels  in  it. 

Similarly,  leftmost  and  rightmost  columns  and  the  bottom  row  are 
tested. 

7)  The  target  is  tested  for  acceptable  size.  The  clusters  must 
be  at  least  3X3  pixels,  and  the  ratio  of  length  of  width  must 

be  between  8:1  and  1:8.  Also,  the  number  of  non-zero  pixels  in  the 
(assumed)  rectangular  boundary  must  be  at  least  75%  of  the  maximum 
possible. 

The  program  to  implement  the  first  four  steps  is  HAM1  as  listed 
in  Appendix  D.  The  results  for  the  infrared  image  in  Figure  20  after 
each  of  the  four  steps  are  shown  in  Figures  21-24.  If  the  last  three 
steps  operated  on  Figure  24,  the  result  would  look  like  Figure  25, 
which  was  altered  manually  for  inclusion  here.  There  are  indeed 
three  targets  as  indicated. 
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Variations  on  the  edging  algorithm  were  tested  but  the  results 
were  the  same  as  with  the  Kirsch  operator.  Instead  of  pursuing  work 
on  the  edging  process,  other  parts  of  the  preprocessing  were  changed 
with  hopes  improving  the  entire  algorithm. 


LOCAL  THRESHOLDING 


With  all  of  the  other  steps  unchanged,  the  conjunctive 
thresholding  was  implemented  using  local  neighborhood  statistics. 

The  threshold  for  each  pixel  is  separately  computed  using  the  mean 
plus  the  standard  deviation  of  the  square  local  neighborhood.  The 
computation  time  increases  significantly,  but  the  process  now  has 
the  ability  to  correct  for  variations  in  the  average  intensity  across 
the  image. 

The  results  for  the  local  thresholding  were  not  very  good  until 
another  step  was  also  changed.  When  using  local  thresholding,  the 
edged  image  must  be  formed  from  the  original  image  instead  of  the 
enhanced  image.  This  is  because  the  enhanced  image  dilates  the  hot 
blobs,  and  when  edged,  the  conjunctive  thresholding  shows  that  the 
thresholded  enhanced  and  edged  images  will  not  coincide  very  well. 

By  edging  the  original  image,  the  edges  of  the  hot  blobs  are  closer 
to  the  center  of  the  cluster.  Then  the  thresholded  edges  coincide 
better  with  the  thresholded  enhanced  image. 

The  program  to  perform  the  new  algorithm  with  local  thresholding 
is  HAM21  as  listed  in  Appendix  D.  The  result  for  this  version  of 
the  algorithm  using  a  17  X  17  neighborhood  for  Figure  20  is  shown 
in  Figure  26.  The  last  three  steps  of  Hamadani's  algorithm  would 
reduce  the  image  to  Figure  27.  The  three  targets  are  again  detected. 


NEIGHBORHOOD  SIZE 


The  amount  of  computation  for  local  thresholds  is  proportional 
to  the  square  of  the  size  of  the  neighborhood.  We  are  thus  motivated 
to  try  to  obtain  good  results  with  smaller  neighborhoods.  17  X  17 
neighborhoods  produce  good  results,  but  7X7  neighborhoods  do  not 
unless  some  changes  are  made  in  the  way  the  thresholds  are  computed. 

Specifically,  the  amount  of  local  standard  deviation  that  is 
added  to  the  local  mean  must  be  modified.  A  coefficient  on  the 
standard  deviation  was  experimented  with,  but  good  results  could  not 
be  produced  without  limiting  the  range  of  the  product  of  the  local 
standard  deviation  and  its  coefficient.  The  limits  must  be  imposed 
because  when  the  small  neighborhood  is  placed  over  a  high  contrast 
edge,  the  standard  deviation  will  be  unappropriately  high.  Thus, 
a  maximum  must  be  found. 

In  addition,  small  neighborhoods  allow  some  false  targets  to 
appear  in  the  results.  This  is  because  in  relatively  edgeless  areas 
the  standard  deviation  is  inappropriately  low.  Thus,  a  minimum 
must  be  found.  When  the  standard  deviation  is  modified  by  a  co¬ 
efficient,  a  minimum,  and  a  maximum,  the  results  for  local  thresholding 
with  7X7  neighborhoods  are  better  than  ever. 

The  program  for  local  thresholding  with  the  modified  standard 
deviation  term  is  HAM31  as  listed  in  Appendix  D.  The  result  of  using 
a  7  X  7  neighborhood  in  HAM31  on  Figure  20  is  shown  in  Figure  28. 

The  last  three  steps  of  Hamadani's  algorithm  would  reduce  the  image 
to  Figure  29. 
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THRESHOLD  RULE 

While  making  the  previously  mentioned  improvements  in  Hamadani's 
algorithm,  the  results  of  thresholding  the  enhanced  image  were 
viewed  separately  from  the  thresholded  edge  image.  This  allows  inspection 

of  what  exactly  the  conjunctive  thresholding  does.  With  local  threshold¬ 
ing,  the  process  was  observed  to  depend  mostly  on  the  thresholding 

of  the  enhanced  image. 

The  edging  step  was  then  omitted  to  test  this  observation  for 
another  shortcut.  The  thresholding  is  now  only  on  the  enhanced  image. 

The  program  to  implement  this  shortcut  is  HAM61  as  listed  in  Appendix 
D.  The  result  of  using  a  7  X  7  neighborhood  in  HAM61  on  Figure  20 
is  shown  in  Figure  30.  The  last  three  steps  of  Hamadani's  algorithm 
would  reduce  the  image  to  Figure  31. 

The  preprocessing  using  Hamadani's  first  four  steps  and  all  of 
the  described  modifications  were  tested  on  eight  different  images. 

The  results  are  included  in  Appendix  E. 


Figure  30  -  Locally  Thresholded  Image  of  I  MAG  1 . 1 R  (Using 
7X7  Neighborhoods  and  No  Edging)  after  Connectivity  Test 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


EDGING  OPERATORS 

Large  masks  for  edging  operators  only  decrease  the  resolution  of 
the  edges  but  may  help  the  overall  process  if  the  images  are  noisy. 

The  best  masks  represent  an  ideal  edge  at  various  orientations.  By 
finding  the  mask  that  produces  the  maximum  convolution  sum,  edging 
determines  the  contrast  and  the  orientation  of  the  edges. 

Mean  deviation  is  very  similar  to  difference,  gradient,  or 
derivative  operators  such  as  Roberts  and  Sobel  operators  for  extract¬ 
ing  edges.  These  are  all  more  susceptible  to  noise  than  the  edge 
mask  approach.  Mean  deviation  is  preferred  to  standard  deviation 
because  it  is  faster  to  compute  and  is  less  affected  by  noise. 

CORRELATION  WITH  EDGES 

Template  matching  can  be  implemented  well  by  correlation  with 
edges.  The  template  is  a  line  trace  drawing  of  the  object  of  interest. 
The  correlation  process  described  in  this  effort  corrects  target 
translation  only.  When  separated  orientation  edge  images  are  used, 
the  correlation  process  is  further  improved.  Some  research  into 
whether  more  orientations  help  the  process  is  suggested.  Also,  an 
orientation  correction  correlation  scheme  might  be  devised  using 
the  separated  orientation  edge  images. 


CLUSTER  RECOGNITION 

The  Hamadani  algorithm  works  well  for  detecting  clusters  in 
the  infrared  images.  Local  thresholding  improves  the  process  by 
compensating  for  changing  average  intensity  across  the  image.  17  X 
17  neighborhoods  prove  to  work  well  for  local  thresholding.  7X7 
neighborhoods  can  be  used  if  the  threshold  computations  are  slightly 
modified.  Less  than  7X7  neighborhoods  could  not  be  implemented 
with  good  results.  When  local  thresholding  is  used,  the  conjunctive 
threshold  rule  can  be  reduced  to  thresholding  the  enhanced  image 
only.  Thus  the  edge  image  can  be  omitted. 

Further  work  on  the  cluster  recognition  algorithm  might 
include  computing  the  thresholds  using  selective  neighborhoods 
that  are  perpendicular  to  the  edges.  Enhancement  might  be 
performed  more  than  once.  The  last  four  steps  of  Hamadani's  algorithm 
are  subject  to  research.  Some  better  connectivity  tests,  thinning 
routines,  and  size  tests  could  be  investigated. 
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SUBROUTINE  I OF ( N, MAIN, FI, F2, F3, MS, SI, S2, S3) 


Written  by  Lt.  Simmons  10  Sep  1981 

Version  2 

This  FORTRAN  3  subroutine  will  reed  from  the  file 
COM.  CM  (FCOM.  CM  in  the  foreground)  the  program  name, 
any  global  switches,  and  up  to  three  local  file 
names  and  corresponding  local  switches. 

Calling  arguments: 

N  is  the  number  of  local  files  and  switches  to  be 
read  from  <F)COM.  CM.  N  must  be  1,  2,  or  3. 

MAIN  is  an  ASCII  array  for  the  main  program  file  name. 

FI,  F2,  and  F3  are  the  three  ASCII  arrays  to  return 
the  local  file  names. 

MS  is  a  two-word  integer  array  that  holds  any  global 
swi  tches. 

SI,  S2,  and  S3  are  two-word  integer  arrays  that 
hold  the  local  switches  corresponding  to  FI  through 
F3  respectively. 


Dimension  the  arrays. 

DIMENSION  MAIN! 7) , MS(2) 

INTEGER  FI < 7) , F2<  7) , F3( 7), SI (2) , S2(2) , S3(2 > 

Check  the  bounds  on  N. 

IF!N.  LT.  1.  OR.  N.  CT.  31ST0P  "N  out  of  bounds  in  IOF .  “ 
Process  the  data  in  <F)COM.  CM 

CALL  GROUND(I)  (Find  out  which  ground  program  is  in 
IF<  I.  EQ.  OJOPEN  0,  "COM.  CM"  iOpen  ch.  0  to  COM.  CM 

IF! I.  EG.  1) OPEN  0, "FCOM.  CM"  (Open  ch.  0  to  FCOM.  CM 

CALL  COMARG! 0, MAIN, MS, IER )  ; Read  from  (F)COM. CM 

IF!IER.NE.  I ) TYPE"  COMARG  error: ”, IER 
WRITE! 10,  1 )MAIN( 1 )  i Type  program  name 

FORMAT!'  Program  ',S13,  'running.  ') 

CALL  COMARC !0, FI, SI, JER)  ; Read  from  !F)COM.  CM 

IF!JER.NE.  I ) TYPE"  COMARG  error  !F1):",JER 

IF!N.  EG.  1)00  TO  2  i Test  N 

CALL  COMARG <0, F2, S2, KER)  (Read  from  !F)COM  CM 

IF (KER.  NE.  1 )TYPE"  COMARG  error  (F2):",KER 

IF!N.  EG.  2) GO  TO  2  (Test  N 

CALL  COMARG !0>  F3, S3, LER )  (Read  from  !F)COM.  CM 

IFILER.  NE.  DTYPE"  COMARG  error  !F3):",LER 

CLOSE  0 

RETURN 


SUBROUTINE  UNPACX(N.  PIXWORD.  PIXELS) 


Written  by  Lt.  Simmons  Version  2 

This  subroutine  mill  unpack  four  4-bit  integers  from  a 
16-bit  integer  word.  The  pixels  in  a  video  file  have  to 
be  unpacked  if  each  pixel  is  to  be  operated  on  separately. 

INTEGER  PIXUORD(N). PIXELS (4. N)  ; Four  pixels  per  word 

DO  1  1*1. N  i  'N'  allows  higher-order 

DO  1  J*l>4  ; arrays  to  be  passed. 

PIXELS< (S-J),  I)*13.  AND.  PIXWORD<I>  ;Pick  off  right  pixel 
PIXWORD< I )*ISHFT (PI XWQRD( I ) . -4 )  iShift  word  4  bits  right 

RETURN  i to  pick  off  next  pixel. 

END 
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Program  PICT  Writtan  by  Lt.  Simmon*  14  Oct  1981 

Fortran  9  Ravisad  by  Lt.  Cromar  12  July  1982 

Ravisad  by  Capt  Wall*  29  Sapt  1983 

This  program  will  convart  vidao  pixels  to  linaprintar  pixels 
and  will  sand  tha  pictura  to  tha  Eclipse's  Printronix  300 
linaprintar.  This  program  prints  tha  complata  296  by  296 
pixal  pictura.  Odd  numbarad  vidao  linas  ara  raprasantad  by  3x3 
pixals.  whila  avan  numbarad  vidao  linas  ara  raprasantad  by  4x3 
pixals.  Run  tima  should  bo  about  two  minutas. 

DIMENSION  1 1 ARRAY (268 ) .  I2ARRAY ( 264 ) • I3ARRAY < 201 > .  I4ARRAY(19S> 
EQUIVALENCE  ( 1 1 ARRAY.  I2ARRAY.  I3ARRAY, I4ARRAY) 

DIMENSION  ILP < 4. 67) .  IFILE<7>.  IREC(64>,  ISAV<4> 

INTEGER  MAIN< 7) >  F2< 7 > .  F3< 7 ) .  MS <  2) »  SI ( 2) >  S2 ( 2 ) < S3< 2) 

CALL  I OF < 1. MAIN. 1FILE.  F2.  F3.  MS. SI. S2. S3) 

C  Sat  up  solid  lina.  spaca.  and  lina  faod/p lot-on  charactar 


IL-177777K  (Solid  lina 

NC=<40100K  i  Spaca 

LF-012K  i  Lina  faed 

LFPC-2412K  iLinx  feed/plot  on 

C  Sat  up  paramatars  for  complata  pictura  display. 

Nl-66  ;  Top  and  bottom  bordar  langth 

N2-296  {Number  of  linas  displayed 

N3»l  ; Location  of  left  bordar 

N4-66  ^Location  of  right  bordar 

N9-67  {Location  of  lina  feed 

N6-1  {Length  of  linas  displayed 


CALL  0PEN(1. IFILE, 1.  IER)  {Open  tha  vidao  file 

CALL  CHECK (IER) 

C  Put  a  bordar  at  tha  top  of  tha  pictura. 

9  DO  7  1-1.  3 

DO  6  J-1.N1 

WRITE  BINARY ( 12 ) IL 

6  CONTINUE 
WRITE  BINARY! 12)LFPC 

7  CONTINUE 

C  Each  lina  of  tha  picture  will  have  a  bordar  on  each  and. 

JTEST  »  -1 

DO  13  JA-1.N2 

JTEST  -  JTEST  *<-!> 

JJ  ■  4 

IF ( JTEST.  GT.  0)JJ»3  {Odd  iteration 

C  Put  a  bordar  down  tha  left  hand  side 

DO  8  K-l.JJ 

ILP  <  K.  N3)-43900K 


{Print  a  lina 
{Terminate  tha  lina 
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CONTINUE 


PAGE 


Put  a  border  down  the  right  hand  side 

DO  9  L-l.JJ 

ILP<L> N4)*40170K 
ILP !L< N3)*LFPC 
9  CONTINUE 

Convert  the  video  picture  pixels  to  ti.PT  pixels 

READ  BINARY! 1 ) IREC  *  Read  one  video  line 

00  12  N-N6.  64 

IWR-BYTE! IRECCN).  2)  iRight  two  pixels 

IWL»BYTE!IR£C!N).  1)  .Left  two  pixels 

IF!JJ.  GT.  3.  3 ) CALL  QUT4X3 ( IWL.  ISAV) 

IF< JJ.  LT.  3.  3 > CALL  0UT3X3(IWL.  ISAV) 

00  11  JB-l.JJ 

BYTE! ILP! JB. N+l )< 1 )*ISAV! JB)  ;  Hove  to  high  byte 

11  CONTINUE 

IF! JJ.  LT.  3.  3) CALL  0UT3X31IWR.  ISAV) 

IFIJJ.  OT.  3.  3) CALL  0UT4X3!IWR.  ISAV) 

DO  12  JC-l.JJ 

BYTE! ILP! JC.  N+l).  2)-ISAV!JC)  ; Store  low  byte 

12  CONTINUE 

L-0 

DO  130  JE-l.JJ 
00  130  JD-1.N3 
L-L+l 

I1ARRAY!L)-ILP! JE. JO) 

130  CONTINUE 

IF!N3.  EQ.  66.  .  AND.  JJ.  EQ.  3)WRITE  BINARY! 12) I4ARRAY 
IF!N3.  EQ.  67.  .  AND.  JJ.  EO.  3)WRITE  BINARY! 12) I3ARRAY 
IFCN3.  EQ.  66.  .  AND.  JJ.  EQ.  4JWRITE  BINARY! 12) I2ARRAY 
IF!N3.  EQ.  67.  .  AND.  JJ.  EQ.  4)WRITE  BINARY! 12) 1 1ARRAY 

13  CONTINUE 

Put  a  border  and  title  at  the  botto*  of  the  picture 

DO  13  JF-1.3 

DO  14  JG«1.N1 

WRITE  BINARY! 12) IL 

14  CONTINUE 
WRITE  BINARY! 12)LFPC 

13  CONTINUE 

WRITE! 12» 16)  j Title  picture 

16  FORMAT! '  '.  13X.  'Signal  Processing  Laboratory.  Air  Force 

♦  Institute  of  Technology.  Wright-Patterson  AFB.  OH  43433<14>'> 


»Print  a  line 
» Terminate  the  line 


; Insert  border 
i line  feed  after 
»  the  picture 


CALL  RESET 
STOP 


i Close  all  channels 
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SUBROUTINE  0UT3X3(VIDPIX. LINEPRINT) 


Written  by  Lt.  J.  H.  Cromer 

This  subroutine  converts  the  video  pixel  velues 
(i.e.  an  integer  value  fro*  0-19)  to  lineprinter 
dot  matrix  form.  A  3x3  array  pattern 
is  formed  by  this  subroutine.  Dot  pattern  texture 
(distribution  of  dots)  and  average  brightness  are 
varied  to  create  16  pseudo-gray  levels.  Odd  numbered 
rows  of  the  picture  created  by  PICTURE  use 
these  3x3  patterns. 

NOTE:  The  six  least  significant  bits  of  each  byte 

sent  to  the  P-300  represent  print  hammer 
switches  (i.e.  a  1  turns  the  hammer  on 
to  print  a  dot.  a  0  leaves  it  off).  Bit 
seven  must  have  a  value  of  1. 

(See  the  Printronix  manual  for  further  discussion) 


>•*****•••****•• *************** ****************************** 

INTEGER  VIDPIX.  LINEPRINT(3> .  RIGHT.  PATTERN (3. 16. 2> 

Note  that  right  and  left  pixel  patterns  are 
not  necessarily  the  same. 

DATA  PATTERN/4*7.  9.  7,  6.  7.  3.  7,  9.  2*2.  7.  2.  9.  2.  9.  2.  7.  2*2. 

«9.  2.  9.  O.  9.  2.  9.  2.  9.  0.  2*2.  0.  2*2.  0.  2.  1.  3*0.  2,  4*0. 

•4*1 70K. 190K. 170K. 160K. 170K, 130K. 120K. 190K. 170K. 4*190K. 

•  120K. 190K.  120K.  170H.  2*120K.  190K. 120K. 190K,  100K.  190K. 

•  140K.  120K.  UOK.  120K.  100K.  190K.  110K.  140K.  UOK.  lOOK.  190K. 

«2*100K. 120K. 7*100K/ 

RI0HT-(VIDPIX.  AND.  19)*1 
LEFT-! ISHFT! VIDPIX.  -4).  AND.  19>>l 
DO  10  1-1.3 

LINEPRINT! I ) -PATTERN! I.  LEFT. 1  >«-PATTERN(  I. RIGHT.  2) 

10  CONTINUE 
RETURN 
END 
C 

C*********  Subroutine  0UT3X3  ************************************ 


oononoonooo 


SUBROUT  X  NE  0UT4X3  <  V I DP  I X , L I NEPR I NT ) 

>•*»••*•****•***•*»*•*•*•***•«***•**«*«** 


Written  by  Lt.  Cromer 

This  subroutine  return*  lineprinter  pixel*  to  the 
celling  program  PICTURE*  which  sends  video 
pixels  (en  integer  from  0-19).  The  pixel  pattern 
returned  is  a  4x3  dot  matrix  array*  to  be  used  for 
the  even  rows  of  pictures  created  by  PICTURE. 

(See  0UT3X3.FR  for  more  explanation). 


C**********  ***********************#*#*#«************************* 

INTEGER  VIDPIX, LINEPR INT( 4 ) .  RIGHT*  PATTERN<4.  16.2) 

DATA  PATTERN/9*7,  9.  7.  7,  6.  7.  7,  3.  3.  6.  7,  9.  9.  6.  3.  2.  9.  2. 

*3*9.  2.  2*  9.  9.  0.  2.  9.  2.  2.  9.  2*  2.  1.  4.  2.  2*  4.  1. 2.  1,  4.  2*  0.  0.  1. 

*4. 0.  0.  2.  10*0. 9*170K.  190M. 2*170K.  160K*  2»170K. 2*130K. 

S160K.  170K* 190K. 120K. 2*170H,  120K.  190K. 120H. 3*190K. 2*120K. 

•  1 90K,  120K.  2*1 90K* 3*120K.  190K* 2*120K.  110K,  140K.  120K.  190K. 
S2*100K.  190K.  120K,  100K.  UOK*  140K.  120K.  2*100K, 

*120K> 2*100K*  140K.  3*100K.  120K.  9*100K/ 

R I OHT “(VIDPIX.  AND.  19M-1 
LEFT-(ISHFT(VIDPIX. -4).  AND.  19)*1 
DO  10  1-1.4 

LINEPR INT<  I ) -PATTERN!  I.  LEFT.  1 )  *PATTERN(  I.  RIGHT.  2) 

10  CONTINUE 
RETURN 
END 
C 

C*********  Subroutine  0UT4X3  ************************************ 


■*** 
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C*a*#a*#a#*e**e*****aa**a*a#****#*#*a*********a**a*#*****a*aa 

Program  TONER  Mrittan  by  Lt.  Jim  Cromar 

This  program  will  convart  individual  pixel  values 
of  an  input  vidao  fila  into  now  valua*  assigned  by 
the  user.  It  can  ba  used  to  adjust  gray-lavals 
(i.a.  histogram  equal i zation ) *  incraasa  contrast, 
display  salactad  pixel  valua*.  craata  "negative  im- 
agasM.  craata  filos  of  a  constant  gray-level.  and  for 
many  othar  purposss.  (NOTE:  Tha  program  EVIOHIST  can 
ba  run  to  ganarata  a  histogram  of  any  VIDEO  fila.  ) 

Exacution  Lina  Format: 

TONER 

Tha  program  will  ask  for  tha  input  and  output 
fila  namas.  and  tha  now  pixel  valuas. 

Load  Lina  Format: 

RLDR  TONER  TIMER  UNPACK  REPACK  CFLIBC 


aaeeeeeaaaaaaaaeaaaeaaaaaeeeeeeaeaeaeaeaeaeeeee** ************* 

INTEOER  NEWVALUE ( 0:  19).  INFILE(7>.  0UTFILE(7> 

I NTEGER  PACKED ( 4096 ) .  UNPACKED  < 1 6384 ) 

ISTART-0  >  start  timar 

I STOP* 1  i stop  timar 


USER  INPUT  OF  VARIABLES  a**************************** 

ACCEPT"Enter  nama  of  vidao  fila  to  ba  tonad  - >  " 

READ (11. 1000 ) INFILE ( 1 ) 

ACCEPT"Entar  nama  of  output  fila - >  " 

R£AD< 11.  1000 >OUTFILE< 1 ) 

TYPE"<19>".  "Enter  now  pixal  vaiues<13>" 

TYPE"  NOTE:  Leading  zeros  are  significant,  i.  a.  enter" 

TYPE"  a  '1'  as  '01',  entar  a  '2'  as  '02'.  ate.  <13>" 

DO  1  J-0. 19  i 16  gray-lavals 

WRITEdO.  2000)  J 
READ  <11. 3000 ) NEWVALUE ( J ) 

IF < NEWVALUE <  J > .  LT.  0.  OR.  NEWVALUE(J).  QT.  13) 
•NEWVALUE (J)™1 3 
1  CONTINUE 

CALL  TIMER! ISTART) 

C 

CMMMM  I/O  FILE  MANAGEMENT  a******************************** 

C 

CALL  OPENU.  INFILE.  1.  IER) 

IFdER.  NE.  1)TYPE" INFILE  OPEN  error  IER 
CALL  DFILW<OUTFILE.  IER) 

IFdER.  NE.  1.  AND.  IER.  NE.  13)TYPE"00TFILE  DFILW 
•  error  IER 

CALL  CFILW<OUTFILE,  3.  64.  IER)  i craata  a  contiguous  fila 
IF< IER.  EQ.  41 >CALL  CFILW(OUTFILE. 2. IER > 


IF! IER.  NE.  1 )TYPE”QUTFILE  CFILW  »rror  #".  IER 
CALL  OPEN (2. OUTFILE. 3. IER) 

IFCIER.  NE.  1) TYPE "OUTFILE  OPEN  irror  IER 
C 

c ********  RE-TONE  THE  PICTURE  ******************************** 

C 

DO  3  J-O.  48.  16 

CALL  RD8LK!  1.  J.  PACKED. 16. IER) 

IF<IER.  NE.  1 )TYPE"RDBLK  #“.J.  “  *rror: ", IER 
CALL  UNPACK<4096. PACKED.  UNPACKED) 

DO  2  1*1.16384  ; do  1/4  of  pictur# 

UNPACKED < I ) -NEWVALUE ( UNPACKED  < I ) ) 

2  CONTINUE 

CALL  REPACK <4096,  UNPACKED.  PACKED) 

CALL  WRBLK<2. J. PACKED, 16. IER) 

IF(IER.  NE.  1)TYPE"WRBLK  #*'.  J>  “  *rror:  IER 

3  CONTINUE 
C 

C  S*nd  ■••tag*  to  CRT  terminal 

C 

CALL  TIMER <ISTOP) 

WRITE! 10.  4000>0UTFILE< 1 ) 

TYPE"<13>". "Hav*  a  nie*  dag!<7><13>" 

C 

C *********  WRITE  NEW  TONE  VALUES  TO  THE  LINEPRINTER  *************** 

C 

WRITE! 12. 3000) 

WRITE! 12.  6000) INFILE! 1 ). OUTFILE! 1 ) 

DO  3  1-0. 13 

WRITE! 12.  7000 ) I . NEUVALUE < I > 

3  CONTINUE 

1000  FORMAT  <S13) 

2000  FORMAT!”  Chang*  old  pix*l  valu*“. 13. "  to  ?”) 

3000  FORMAT! 12) 

4000  FORMAT!"  Th*  ton*d  pietur*  i*  in  th*  fil*  - >  “,S13) 

3000  FORMAT! //////26X.  "  RESULTS  OF  TONER<10>"/ 

•26X.  " - "////) 

6000  FORMAT ! 10X, "  Input  fil*  - >  ", S13. /10X, "  Output  fil* 

• - >  ",  313. //20X.  "OLD  PIXEL".  10X. "NEW  PIXEL". /20X, 

«" - ",  10X.  " - "/) 

7000  F0RMAT!23X.  12,  17X.  12) 

CALL  RESET 
STOP 
END 
C 

C ***********  Program  TONER  **•••**•***«••**•******«*«**««*•**•***«**« 
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SUBROUTINE  TIMER (I ) 


■•»**■**#****#*■****#*#***##**#***•*•***■***•**•*•*#■*■****•*•*•**#***•*•*** 

Subroutine  TIMER  Written  by  Lt.  Jim  Cromer 

Fortran  S 

Thie  subroutine  is  used  to  time  the  reel-time  execution 
time  of  the  celling  program.  If  the  parameter  passed.  I. 
is  equal  to  0.  the  timer  is  unconditionally  started. 

If  I  is  not  equal  to  O.  the  timer  is  unconditionally 
stopped,  and  the  total  run  time  is  typed  on  the  console 
CRT. 

Execution  Line  Format 

CALL  TIMERd)  ilFd.EQ.  0).  start  timing 
;IFd.NE.  O).  stop  timing 

Ct#M*M*M**«*M««4**«»********t**M«*««M4«**««*M*t<«*»**«******4 

COMMON  /ITIME/  1H1,  IM1.  IS1 
IFd.  NE.  0)00  TO  100 

CALL  FQTIMEdHl,  IH1.  IS1>  »  get  starting  time 

WRITEdO.  1000)  IH1.  IM1.  IS1 

lOOO  FORMAT ( /  /"  START  TIME - 14.  -.  13.  “.  13) 

RETURN 

too  CALL  F0T1ME11H2. ina. XS2)  > get  stopping  time 

WRITE < 10. 2000) IH2.  IM2.  IS2 

2000  FORMAT < //“  STOP  TIME - >".  14.  ",  13,  ".  13) 

IT0TAL-3600*< IH2-IM1 )+60*< IM2-IM1 ) +IS2-IS1 
HOURS- 1 NT  < IT0TAL/3600 ) 

TRQN— < IT0TAL-3600*HQURS>  i intermediate  variable 

M I NS— I NT  <  TR0N/60 ) 

I  SECS— MOO  <  TRW.  60) 

WRITEdO.  3000)HOURS.  MINS.  ISECS 

3000  FORMAT*//"  TOTAL  TIME - >", 14.  ",  13,  ".  13) 

RETURN 
END 
C 

Ceeae**********  Subroutine  TIMER 
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Program  NMOVE  Written  by  Lt.  Jim  Cromor 

Fortran  IV  16  Aug  1982 

This  program  will  plaea  tha  vidto  information  in  tha  win¬ 
dow  givan  for  tha  tamplata  (inputfilal)  insida  of  tha 
window  givan  for  tha  background  ( inputf ila2)»  and  writa 
tha  combinad  pictura  to  tha  outputfila.  Tha  window  may  ba 
placad  anywhara  within  tha  background>  and  may  ba  takan 
from  anywhara  within  tha  tamplata.  Window  width,  langth. 
and  position  ara  input  by  tha  usar. 

Eiacution  lina  format:  (on  tha  NOVA  only) 

NMOVE  (run  EMOVE  on  tha  ECLIPSE) 

Loadar  command  lina  format  (NOVA  only): 

RLDR  NMOVE  TEST  BLOCK  CHANGE  XWRBLK  XRDBLK 
UNPACK  REPACK  FORT.  LB 


INTEGER  IPAR (2) >  INFILE1 (7). INFILE2(7) .  0UTFIL£(7> . 

SCB1. C BLOCKS.  CCOL.  CLS.  CSTOP.  CTOP.  CLB.  CLEFT.  TTOP.  TB1.  TBLOCKS.  CH3. 
STCQL.  TLS. TSTOP.  TLB.  TLEFT,  COMB ( 1024) . TEMP ( 1024) , BACK ( 1024). WIDTH. TB 
COMMON/LI ST 1 /COMB.  TEMP.  CLS.  TLS 
COMMON/L I ST2/LENGTH.  WIDTH 
EQUIVALENCE  (COMB.  BACK) 


**  I/O  FILE  MANAGEMENT  a********************************** 

TYPE"<19>". "Program  NMOVE  is  to  ba  run  on  tha  NOVA  only!" 

99  TYPE"<19>".  "**********•*♦******#•-*******•»**•»*■***-*■**♦•»♦***<  13>“ 

ACCEPT"  Entar  tamplata  fila  nama:  ” 

READ (11. 1000) INFILE1 ( 1 ) 

ACCEPT"<1S>".  "  Entar  background  fila  nama:  " 

READ (11. 1000) INFILE2( 1 ) 

DO  999  J-1.7 

999  0UTFILE(J)-INFILE2(J> 

ACCEPT"<19>".  "  Entar  combinad  output  fila  nama:  " 

READ (11. lOOO)OUTFILE(l) 
lOOO  FORMAT ( S13) 

CALL  0PEN(1.  INFILEl.  2.  IER) 

IFdER.  NE.  1)TYPE"  Channal  1  OPEN  arror:  ".  IER 
CALL  OPEN ( 2.  INFILE2.  2.  IER) 

IFdER.  NE.  1) TYPE"  Channal  2  OPEN  arror:  ".  IER 
CH3-2 
I COUNT "O 

DO  1002  J>1.7  ichack  for  BACKGROUND-COMBINED 

1002  IF(OUTFILE(J).  EQ.  INFILE2( J) ) ICOUNT-ICOUNT+1 

IFdCOUNT.  EQ.  7)00  TO  1 
CH3—3 

CALL  DFILW(OUTFILE.  IER) 

IFdER.  NE.  1.  AND.  IER.  NE.  13)TYPE"0UTFILE  DFILW  arror:  ",  IER 
CALL  CFILW(OUTFILE.  2.  IER) 

IFdER.  NE.  DTYPE"  CFILW  arror:  ".  IER 
CALL  QPEN(CH3< OUTFILE.  2.  IER) 

IFdER.  NE.  DTYPE"  Channal  3  OPEN  arror:  ",  IER 


A-  12 


PACE  2 


C*******  ENTER  WINDOW  PARAMETERS  ******************************* 

C 

1  ACCEPT"<15>"»  "  Enter  top  row  of  template  window  ( 1-256) :"» TTOP 

ACCEPT"  Enter  left  column  of  template  window  ( 1-256) :". TLEFT 
ACCEPT"  Enter  width  of  window  ( 1-256) WIDTH 
ACCEPT"  Enter  length  of  window  ( 1-256) LENGTH 
C 

C  The  call*  to  TEST  check  to  see  if  the  input  parameters  are 

C  legal,  and  modifies  them  if  necessary: 

C  0  <  TOP  <  257.  (TOP  ♦  LENOTH)  <  258 

C  0  <  LEFT  <  257.  (LEFT  *  WIDTH)  <  258 


CALL  TEST (TTOP.  TLEFT) 

ACCEPT"<15>". "  Enter  top  row  of  background  window  ( 1-256) CTOP 
ACCEPT"  Enter  left  column  of  background  window  ( 1-256) CLEFT 
CALL  TEST (CTOP.  CLEFT) 

CALL  BLOCK (T8L0CKS.  TB1.  TLS.  TCOL,  TTOP. TLEFT) 

CALL  BLOCK (CBLOCKS.  CB1.  CLS.  CCOL.  CTOP.  CLEFT) 

C 

C  Determine  column  number  of  the  last  video  row  (0-3) 

C 

J1 -MOD (LENGTH, 4) 

TSTOP— MOD ( ( TCOL* J 1 > .  4 ) - 1 
CSTOP— MOD ( ( CCOL+J 1 ) .  4 ) - 1 
IF(CSTOP. EQ.  -1 ) CSTOP— 3 
IF (TSTOP.  EQ.  -1 ) TSTOP— 3 
C 

C  Determine  the  last  significant  block  of  window 

C 

CLB— CB 1 ♦CBLOCKS— 1 
TLB— TB 1 +TBLOCKS- 1 
C 

C  User  check  of  window  parameters 

C 

TYPE"<15>“.  "WIDTH-".  WIDTH.  "  LENGTH-". LENGTH 

TYPE" TEMPLATE  TOP  ROW-", TTOP, "  BACKGROUND  TOP  ROW-",  CTOP 
TYPE " TEMPLATE  LEFT  COLUMN-".  TLS.  "  BACKGROUND  LEFT  COLUMN 
4— ".  CLS.  "<15>" 

ACCEPT'Enter  1  to  see  expanded  set  of  variables,  ang 
•  other  integer  to  continue:  ".  I 
IF(  I.  NE.  1)00  TO  5 


TYPE"************************************************" 

TYPE"  PARAMETER  TEMPLATE  BACKGROUND" 

TYPE"  _  _  _ " 

TYPE"<15>"» "  TOP  ROW  ".TTOP, CTOP 
TYPE"<13>", *  START  COL  TCOL,  CCOL 
TYPE"<15>". "  STOP  COL  •  ".  TSTOP.  CSTOP 
TYPE"<15>"» "  FIRST  BLOCK". TB1. CB1 
TYPE"<15>",  "  LAST  BLOCK  ".TLB, CLB 
TYPE"<15>". "  •  OF  BLOCKS".  TBLOCKS,  CBLOCKS 
TYPE"<15>"» "  LEFT  COL  ",TLS. CLS 
TYPE"<15>". "  WIDTH-  “.WIDTH 
TYPE"<15>". "  LENGTH-".  LENGTH 

TYPE " < l 5> " . " *********************************************** " 
ACCEPT"  Enter  1  to  try  another  set,  any  other  integer 
t  to  continue:  ". I 
IF( I. EQ.  1)C0  TO  l 


A-13 


UOO  OOUU  UUUOOUUOOOuu  uuo 


PACE 


3 


C**##**  Create  th*  combined  picture  **************************** 

C 

5  I COUNT-O 

IF(CH3.  EG.  2) CO  TO  20  ilf  combined  picture  file 

> is  the  same  as  the  back- 
j ground  picture  file,  then  no 
i  need  to  write  to  itself 
Write  background  only  blocks  (before  window) 
to  the  combined  picture  file. 

JMAX-CBl-l 

IF( JMAX.  LT.  0)00  TO  20 
DO  10  J-O,  JMAX 

CALL  RDBLK ( 2< J. BACK. 1. IER) 

IFUER  .  NE.  DTYPE"  2RDBLK".  J.  "  error:  IER 
CALL  URBLK(CH3> J.  COMB.  1.  IER) 

IF(I£R.  NE.  1>TYPE"  WRBLK". J. "  error:  ",  IER 
IC0UNT-1C0UNT+1 

10  CONTINUE 

20  TYPE"  Background  before  window  completed.  “ 

TYPE"  #  Blocks  written: ".  ICOUNT 


Overlay  template  window  onto  background  . 

CALL  XRDBLK< 1.  TB1.  TEMP.  1,  IER) 

IF< IER.  NE.  1) TYPE" 1 RDBLK  TB1, “  error: ", IER 

CALL  XRDBLK < 2.  CB1.  BACK.  1.  IER) 

IF (IER.  NE.  1 )TYPE"2RDBLK  #".CB1,“  error: ". IER 
Nl-TCOL  i  4-MAX <Nl,  N2)  gives  the  number  of  rows 

N2-CC0L  s to  change  before  the  next  RDBLK 

IF(TCOL.  OT.  CCOL)CO  TO  lOO 


There  are  four  columns  in  the  packed  video  array  (64x4). 
designated  O,  1.  2,  and  3.  If  the  template  starting 
column  number  is  less  than  or  equal  to  the  background  (combined) 
starting  column  number,  then  the  background  block  will  be  ‘used 
up"  before  the  template  block.  When  the  background  block  is 
finished,  a  WRBLK  is  done,  and  the  next  background  block  is  read. 
When  the  template  block  is  finished,  the  next  template  block  is 
read,  but  no  WRBLK  needs  to  be  performed.  Note  that  the  back¬ 
ground  and  combined  files  are  always  at  the  same  block  number. 

CALL  CHANCE (N2.  N2.  N1 ) 

CALL  XWRBLK( CH3.CB1,  COMB.  1.  IER) 

IF( IER. NE. 1) TYPE"  WRBLK  »".CBi."  error: ", IER 

Write  the  template  window  into  the  background 

TB-TB1+1 
IMIN-CBl+1 
I COUNT- 1 
IMAX-CLB-1 

IFtIMIN.  CT.  I MAX ) 00  TO  60 
DO  90  I-IM1N.  IMAX 

CALL  XRDBLK (2,  I.  BACK, 1, IER) 

IF( IER.  NE.  1) TYPE*  2RDBLK  error:  ", IER 

CALL  CHANGE (Nl. N2. N1 > 

CALL  XROBLKU.TB.  TEMP,  1.  IER) 

IF(  IER.  NE.  DTYPE"  1  RDBLK  #".  TB.  •  error:  ",  IER 
CALL  CHANOE ( N2,  N2, Nl ) 
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CALL  XWRBLK < CH3.  I, COMB.  1,  IER ) 

IF <  IER.  NE.  1  )TYPE’*  WRBLK  I,  "  error:  ",  IER 
I COUNT* I COUNT fi 

30  TB*TB+1 

60  TYPE"  TCOL.  LT.  CCOL — Window  portion  complete." 

TYPE"  •  block*  written: ICOUNT 
GO  TO  230 
C 

C  In  this  case  the  template  starting  column  number 

C  is  greater  than  the  background  starting  column  number. 

C  The  template  block  must  be  "finished"  first. 

C 

100  CALL  CHANCE <N1.  N2.  N1 )  ; finish  TEMP  block 

TB-TB1+1 
IMAX-CLB-1 
I COUNT-O 

IF  ( CB1.  CT.  I  MAX  )  GO  TO  223 
00  200  I*CB1. IMAX 

CALL  XRD8LK< 1. TB. TEMP.  1. IER) 

IF< IER. NE.  1 )TYPE"  1RDBLK  •  ". TB.  "  error:  IER 
CALL  CHAN0E(N2, N2, Nl>  {finish  BACK  block 
CALL  XWRBLK(CH3,  I. COMB.  1.  IER) 

IF <  IER.  NE.  1 )  TYPE"  WRBLK  •  ’*,  I.  "  error:  IER 

ICOUNT*ICOUNTrl 

IBLK* I +1 

CALL  XRDBLK ( 2.  IBLK. BACK.  1,  IER) 

IF< IER. NE.  1 )TYP£“  2R0BLK  •  “.  IBLK. «  error:  “.  IER 
CALL  CHANGE(N1.  N2. N1 )  {finish  TEMP  block 

200  TB-TB+1 

223  TYPE"  TCOL.  OT.  CCOL—Window  portion  complete." 

TYPE"  •  blocks  written: *. ICOUNT 
C 
C 

C - 

c 

C  If  the  combined  <background>  stopping  column  number 

C  is  greater  than  the  template  stop  column  number,  then  the 

C  second  last  template  block  <I2LB*TLB-1>  must  be  read  (i.e 

C  there  are  more  video  rows  to  be  changed  in  the  last  back- 

C  ground  block  than  there  are  available  in  the  last  template 

C  block  to  change  them  to).  If  TSTOP  is  greater  than  or  egi<al 

C  to  CSTOP.  then  there  are  sufficient  video  rows  available  in 

C  the  last  template  block  to  complete  the  last 

C  background  block  to  be  changed. 

C 

230  IFiCSTOP.  OT.  TSTOP) CO  TO  400 
Ml-TSTOP -CSTOP 

CALL  XRDBLK ( 1. TLB.  TEMP.  1, IER) 

IF< IER.  NE.  1 )TYPE"  1RDBLK  •".TLB. ■  error: ", IER 
CALL  XRDBLK <2. CLB. BACK. 1, IER) 

IF< IER.  NE.  1) TYPE"  2RDBLK  *".CLB. "  error: ", IER 

N1-3-CST0P 

N2-0 

CALL  CHANCEfNl,  N2.  Ml)  {finish  BACK  block  to  CSTOP 
CALL  XWRBLK < CHS, CLB. COMB. 1. IER) 

IF< IER.  NE.  1 )TYPE"  WRBLK  •“.CLB."  error: ", IER 

TYPE"  CSTOP,  LT.  T8T0P-—  Last  block  of  window  complete.  " 

00  TO  300 
C 

C  Complete  the  last  block  of  the  window 

C  NOTE:  CSTOP  is  greater  than  TSTOP.  Therefore  finish 

C  TEMP  before  BACK. 

C 

400  Ml “CSTOP -TSTOP 
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PACE 


3 


CALL  XRDBLK(2.  CLB.  BACK,  1.  IER) 

IF< IER. NE.  1 )TYPE"  2R0BLK  #",CLB,"  error:  ",  IER 
I2LB-TLB-1 

CALL  XRDBLKU,  I2LB.  TEMP,  1,  IER> 

IF< IER.  NE.  1 ) TYPE"  1RDBLK  I2LB,  “  error:  ", IER 

N1-4-M1 

N2-0 

CALL  CHANCE(N1,N2, Nl)  j finish  TEMP  block 

CALL  XRDBLKU.  TLB,  TEMP,  1.  IER) 

IF< IER.  NE.  1 >TYPE"  1RDBLK  TLB,  “  error: ”, IER 
M1-3-TST0P 

CALL  CHANCE ( Ml, N2, Nl )  {finish  BACK  block  to  CSTOP 

CALL  XWRBLKICH3. CLB, COMB. 1.  IER) 

JFUER.NE.  1)TYPE"  WRBLK  #",CLB. "  error:  ",  IER 

TYPE"  CSTDP.  OT.  TSTOP-— Last  block  of  window  complete.  " 

I COUNT- 1 

Finish  the  combined  file  (background  only  portion) 

JHIN— CLB+1 

IF(JM1N.  CT.  63) GO  TO  601  i  if  finished  STOP 

IF(CH3.  EQ.  2) GO  TO  601  iif  COMB INED-BACKCROUND,  STOP 

DO  600  J-JMIN. 63 

CALL  RDBLK<2, J, BACK, 1, IER) 

IF(IER.  NE.  1)TYPE"  2RDBLK  J, H  error:  ",  IER 
CALL  URBLK<CH3, J,  COMB, 1.  IER) 

IF( IER.  NE.  1) TYPE"  WRBLK  J. "  error:  ", IER 
ICOUNT— I COUNT+1 

CONTINUE 

TYPE"  Finished  background  only  portion.  " 

TYPE"  #  blocks  written:  ",  ICOUNT 

TYPE"<13>“. "<7>", "<13>". “  Program  NMOVE  execution  completed.  <T>" 
WR I TE  < 10. 2000)0UTFILE( 1 ) 

FORMAT < "  The  combined  picture  is  in  the  file  — >  ",S13) 


******  Present  Option  Menu  ********************************* 


GO  TO  2010 

!  TYPE"<13>",  "Input  error  Try  again." 

I  TYPE " < 1 3> " . " ******************************************  » 

TYPE"<13>". "What  n*xt?<13>", "Her*  are  the  options:" 

TYPE"<13>"» "<1 1>1  -  Try  another  set  of  window  values" 

TYPE"<13>". "<1 1>2  -  Start  over  with  new  input  pictures” 

TYPE"<13>".  "<1 1>3  -  Display  combined  picture  on  the  video  monitor" 
TYPE"<13>".  "<1 1>3  -  Save  combined  picture  and  ST0P<13>" 

ACCEPT" <ll>En ter  option  - >  ", IOPT 

IF ( IOPT.  EQ.  1  )00  TO  1 
CALL  RESET 

IF ( IOPT.  EQ.  2>G0  TO  99 
IF ( IOPT.  EQ.  3) STOP 

TYPE"<13>". "Check  monitor  -  -  Press  green  CHOPS  control 
•button  to  continue. " 

I DC NT— 4 
IPAR<  1  >-99«?9 
IPAR<2)-0 

WRITE! 10. 3000  >  OUTF I LE ( 1 ) 

3000  FORMAT! "0".  "Picture  being  displayed  - >  ",S13> 

CALL  CHANNELtO, 0. 3. 0, 0. "A". 0. 0. 0. IE. IS)  ;call  abort 

CALL  CHANNEL (3. 1.2. 1. IDCNT. OUTFILE, 64, 0. IPAR, I ERR, ISYS) 

CALL  ERCHK! I ERR,  1,  IDCNT,  1. ISYS) 

TYPE"<13>".  "CHANNEL  currently  not  loaded." 

TYPE”Use  VIDEO  to  display  combined  pictures.  <13>" 

CALL  OPEN! 1,  INFILE1. 2.  IER)  ; r*-OPEN  channels 

IF! IER.  NE.  1 ) TYPE"CH1  RE-OPEN  ERROR:  ",  IER 
CALL  OPEN! 2, INFILE2, 2. IER) 

IF! IER.  NE.  1 ) TYPE"CH2  RE-OPEN  ERROR:  ",  IER 
IF! CH3.  EQ.  3 ) CALL  OPEN! 3.  OUTFILE. 2.  IER) 

IF! IER.  NE.  1 ) TYPE"CH3  RE-OPEN  ERROR:  ",  IER 
00  TO  2010 
END 
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SUBROUTINE  TEST (TOP. LEFT) 


Subroutine  TEST  checks  to  see  if  the  input  parameters 
to  program  NMOVE  are  legal,  and  modifies  them  if 
necessary.  (It  is  also  called  by  DISTANCE.  ) 


INTEGER  TOP, WIDTH 

COMMON/LIST2/LENGTH,  WIDTH 

IF  (LEFT.  LT.  1.  OR.  LEFT.  CT.  256>LEFT-i 

MAXWIDTH-237-LEFT  i picture  has  236  columns 

IF( WIDTH.  GT,  MAXWIDTH. OR.  WIDTH.  LT.  1 >WIDTH-MAXWIDTH 
IF ( TOP.  LT.  1.  OR.  TOP.  GT.  236>TOP»i 

MAXLENCTH-237-T0P  (picture  has  236  rows 

IF  (LENGTH.  GT.  MAXLENCTH.  OR.  LENGTH.  LT.  1)  LENCTH-riAXLENCTH 
RETURN 
END 


••««••*»*  Subroutine  TEST 


■*** 
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SUBROUTINE  BLOCK (NUMBLOCKS.  BL0CK1,  LEFTSIDE.  COLUMN.  TOP.  LEFT) 


Subroutine  BLOCK  determines  the  total  number  of  blocks  to 
be  read  into  the  window,  the  first  block  to  be  read, 
and  the  first  video  row  "column"  number.  This  subroutine 
is  called  by  NMOVE  and  DISTANCE. 


"»«•«*•«•**»•«•*******•**«*•#••***«* 

INTEGER  BLOCK 1. COLUMN. TOP.  REMAINDER. WIDTH 
COMMON/LI ST2/LENCTH,  WIDTH 

BLOCKl-INT( (TOP-1 )/4.  0)  ;  4  rows  per  block 

COLUMN-MOD ( (TOP-1).  4) 

LEFTSIDE-LEFT 
REMAINDER-MOD (LENGTH. 4) 

K 1 -LENGTH+3 
NUMBLOCKS— INT (Kl/4.  0) 

IF < REMAINDER .  Ed.  2.  AND.  COLUMN.  GT.  2 )NUMBL0CKS«NUMBL0CKS+1 
IF ( REMAINDER.  Ed.  3.  AND.  COLUMN.  GT.  1 > NUMBLOCKS-NUMBLOCKS+1 
IF (REMAINDER.  Ed.  0.  AND.  COLUMN.  GT.  0 ) NUMBLOCKS-NUMBLOCKS* 1 
IF (NUMBLOCKS.  GT.  1 ) RETURN 

TYPE "WARNING:  *  Blocks  to  be  read  -".NUMBLOCKS 

PAUSE 
RETURN 
END 
C 

C ************  Subroutine  BLOCK  ************************************* 
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SUBROUTINE  CHANCE ( JMIN, CSTART. TSTART) 


Written  bg  Lt.  Jim  Cromer 

Subroutine  CHANCE  changes  the  corresponding  background 
(i.e.  the  combined  picture)  pixels  to  template  pixels; 
it  is  called  bg  the  program  NMGVE. 


INTEGER  COMB < 1024), TEMP (1024 >. CLS,  TLS.  CSTART. TSTART. WIDTH 
COMMON  /LIST2/  LENOTH. WIDTH 
COMMON  /LIST1/  COMB. TEMP. CLS. TLS 
C 

DO  2  J-UMIN.  3 

K-TSTART*256+TLS  l Set  left  side  of  input ( temp  late ) 

M-CSTART*236+CLS  .and  output  (combined)  windows. 

KMAX-K+WIDTH-1  » Change  values  over  the  width  of  window 
DO  1  L-K,  KMAX 

COMB <M) -TEMP <L) 

1  M-M+l 
TSTART-TSTART+1 

2  CSTART— CSTART+1 

1F(CSTART.  EC.  4)CSTART-0  ‘ reset  row-pointer  if  necessarg 

IF (TSTART.  EC.  4)TSTART»0 
RETURN 
END 
C 

C *********  Subroutine  CHANCE  ***************************************** 
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SUBROUTINE  XRDBLK (CH.  J,  FILE.  I.  IER  > 


by  Lt.  Jim  Cromer 

Subroutine  XRDBLK  perform*  e  RDBLK  to  the  designated 
channel<  read*  a  packed  video  file  block!  and 
returns  an  unpacked  array. 

#****##*******#****###**#**#********♦*#*#********##**#**#*#*##*#*♦♦' 

INTEGER  CH. FIL£< 1024),  VIDE0<23&> 

K«236*I 

IF<  J.  GE.  0.  AND.  J.  LE.  63)  CO  TO  1 

TYPE"ERRQR:  <7>BL0CK  POINTER  OUT  OF  BOUNDS  IN  XRDBLK H 

TYPE"  J*“,J 

STOP 

1  IF<  I.  EG.  1)60  TO  2 

TYPE "ERROR  IN  <7>XRDBLK" 

TYPE"  •  Blacks  to  be  read  -".I 
STOP 

2  CALL  RDBLK (CH.  J.  VIDEO.  I.  IER) 

DO  3  L-l.K 

DO  3  H-1, 4 

ICOUNT-3— M+<L-1)*4 

FILE( IC0UNT)»13.  AND.  VIDEO(L) 

VIDEO<L)«ISHFT(VIDEO<L).  -4) 

3  CONTINUE 
RETURN 


Subroutine  XRDBLK  ****■*♦♦*•*•*•*■*»♦*******•*•*■**■*•**•*■**-*■■**♦ 
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SUBROUTINE  XWRBLK(CH, J.  FILE,  I, IER) 

Qiven  an  unpacked  video  file.  Subroutine  XUR8LK 
will  pack  the  array  and  write  it  to  the  deeignated 
unit  number. 

INTEGER  CH. FILE< 1024).  VIDEO* 2S6) 

K«236*I 

CALL  REPACK(K, FILE.  VIDEO) 

CALL  WRBUMCH,  J,  VIDEO.  I.  IER) 

RETURN 

END 


'  «'■  -'  /•  W  - 
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C  PROGRAM  NCOMB  -  DG  FORTRAN  9  -  BY  CAPT  ROBERT  WELLS.  DEC  1983 

THIS  PROGRAM  COMBINES  UP  TO  FOUR  VIDEO  IMAGES  (PACKED  VIDEO 
FORMAT).  EACH  PIXEL  OF  THE  INPUT  IMAGES  IS  TESTED  TO  SEE  IF 
IT  IS  NON-ZERO.  THE  VALUE  OF  THE  FIRST  NON-ZERO  PIXEL  OF  THE 
INPUT  IMAOES  IS  ASSIGNED  TO  THE  OUTPUT  IMAGE  AT  THAT  PIXEL 
LOCATION.  IF  ALL  INPUT  IMAGES  ARE  ZERO.  THE  OUTPUT  IS  ZERO. 
THE  OUTPUT  IS  STORED  IN  PACKED  VIDEO  FORMAT. 

INTEGER  INFNM1 (7).  1NFNM2<7>.  INFNM3! 7 ) »  I NFNM4 < 7 ) .  OUTFNM ( 7 ) 
INTEGER  INI (2048).  IN2(2048 > . IN3 ( 2048 ) .  IN4( 2048) .  OUT (2048 ) 

I NTEGER  TEMP (912).  FLCNT 

1  FORMAT (S13) 

2  ACCEPT-'ENTER  NUMBER  OF  FILES  TO  BE  COMBINED  ->".  FLCNT 
IF(  (FLCNT.  LT.  2).  OR.  (FLCNT.  GT.  4)  >00  TO  2 

GO  TO  (2,4.3)  FLCNT 

ACCEPT "ENTER  INPUT  FILENAME  ->H 
READ(11, 1)1 NFNM4 ( 1 ) 

CALL  OPENO.  INFNM4.  1.  IER) 

CALL  CHECK (IER) 

3  ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ(11.  1 ) INFNM3( 1 ) 

CALL  OPEN <4.  INFNM3.  1.  IER) 

CALL  CHECK (IER) 

4  ACCEPT "ENTER  INPUT  FILENAME  ->" 

READ (11,  l)INFNM2d) 

CALL  OPENO.  INFNM2.  1,  IER) 

CALL  CHECK! IER) 

ACCEPT “ENTER  INPUT  FILENAME  ->" 

READ! 11.  1 ) INFNM1 ( 1 ) 

CALL  OPEN ( 2,  INFNM1,  1.  IER) 

CALL  CHECK! IER) 

ACCEPT "ENTER  OUTPUT  FILENAME  ->" 

READ (11,  1 ) OUTFNM! 1 ) 

CALL  OPENd.  OUTFNM.  3.  IER  > 

CALL  CHECK! IER) 

DO  14  1-0. 31  i  COUNT  FOR  2K  BUFFER  LOADS 

00  TO  (2,6.9)  FLCNT 

CALL  RD8LKO,  2*1,  TEMP.  2,  IER) 

CALL  CHECK! IER) 

CALL  UNPACK (912,  TEW.  IN4) 

9  CALL  RDBLK (4, 2*1.  TEMP,  2.  IER) 

CALL  CHECK! IER) 

CALL  UNPACK (912.  TEMP,  IN3) 

6  CALL  RDBLK (3.  2*1.  TEMP.  2.  IER) 

CALL  CHECK! IER) 

CALL  UNPACK! 912.  TEMP,  I N2> 
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PACE  2 

CALL  RDBLK<2,  2*1, TEMP, 2'  IER  ) 

CALL  CHECK < IER ) 

CALL  UNPACK! 312.  TEMP.  INI) 

00  13  J-1,2048  »  COMBINING  LOOP 

GO  TO  (2.8.7)  FLCNT 

IF < IN4( J) .  ME.  0)00  TO  9 

7  IF(IN3(J).  NE.  0)G0  TO  10 

8  IF( IN2< J).  NE.  0)G0  TO  11 

IF  <  INK  J),  NE.  0 )  GO  TO  12 

OUT ( J ) “O  i  ZERO  IF  ALL  IN-FILES  ARE  ZERO 
GO  TO  13 

9  OUT ( J)*IN4! J) 

GO  TO  13 

10  OUT ( J) *IN3( J) 

GO  TO  13 

11  0UT!J)-IN2!J) 

GO  TO  13 

12  OUT  <  J)*IN1 < J) 

13  CONTINUE 

CALL  REPACK <312.  OUT,  TEMP) 

CALL  UR8LK( 1.  2*1.  TEMP.  2,  IER) 

CALL  CHECK! IER) 

14  CONTINUE 
CM  L  OCCCT 

STOP " <7><7><7><7>NC  OMB  « 

END 


Notts  on  program  INVERSE 


Use:  Two-d imansional  invars*  discrata  fourier  transform 

Execution  format:  From  CLI>  command  lina  is: 

INVERSE  inputfila/I  C output f i 1  a/03  #/N 
Paramatar  exp  lanations: 

Inputfilt:  The  input  fila  must  b*  in  binary  format,  since  a 

RDBLH  is  usad  to  raad  data  from  tha  file.  Tha  file  it 
assumed  to  contain  complex  data,  must  be  square,  and  must 
have  a  minimum  length  of  64.  If  the  N  by  N  array  has  N  <  64. 
the  columns  should  b*  augmented  with  zeros  so  that  each  column 
can  be  stored  as  on*  block  of  data.  Thus.  Tor  N  <  64.  the 
size  of  the  input  array  must  b*  64  rout  by  N  columns  of  complex 
DFT  samples.  The  origin  of  tha  frequency  axes  must  be  at 
(N/2  1.  N/2  +1). 

□utputfila:  The  output  file  specification  is  optional.  If 

no  output  file  is  specified,  the  inverse-transformed  data  is 
written  back  into  the  input  file,  overwriting  it.  The  output 
data  is  complex  and  binary.  Note:  this  program  runs  MUCH  faster 

if  Outputfil*  is  NOT  specified. 

tt:  An  integer  must  be  specified  to  indicate  the  size  of  the 

input  array  (N  by  N).  which  is  the  size  of  the  DFT  computed.  N 
must  be  a  power  of  two.  N  <  512. 

This  program  can  be  called  from  a  DC  fortran  program  since  its 
last  line  is  a  call  to  FBACK.  For  more  information,  see  the 
DC  Fortran  IV  User's  Manual,  chapter  II-5. 
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Notes  on  program  DIRECT 

Use:  Two-dimensional  discrete  fourier  transform 

Execution  format:  From  CLI<  command  line  is: 

DIRECT  inputf ile/1  #/N  Coutputf ile/03 
Parameter  exp lanations: 

Inputfile:  The  input  file  must  be  in  binary  format  since  a 

RDBLK  is  used  to  read  data  from  the  input  file.  The  input 
file  is  assumed  to  contain  complex  data>  with  the  imaginary 
part  zero.  The  input  array  is  assumed  to  be  square.  with 
a  minimum  row  or  column  length  of  64.  If  an  array  of  size 
N  <  64  is  to  be  input<  each  column  should  be  stored  as  one 
block  of  data  <i.  e.  64  complex  points)  even  though  only 
the  first  N  complex  numbers  in  the  block  correspond  to  actual 
data  points.  Only  N  rows  need  be  stored.  Thus,  for  N  <  64. 
the  size  of  the  input  array  must  be  64  by  N  columns. 

tt:  An  integer  must  be  specified  to  indicate  the  size  of  the 

input  array  <N  by  N).  which  is  the  size  of  the  DFT 
computed.  N  must  be  a  power  of  two.  N  <  512. 

□utputfile:  The  output  file  specification  is  optional.  If 

no  output  file  is  specified,  the  transformed  data  will  be 
written  back  into  the  input  file,  destroying  the  input  data. 

The  output  file  is  complex,  and  is  in  binary  format.  The  origin 
of  the  DFT  samples  will  be  at  location  (N/2.  N/2),  and  the  rows 
and  columns  will  be  transposed.  Note:  this  program  runs  MUCH 
faster  if  Outputfile  is  NOT  specified. 

This  program  can  be  chained  from  a  DO  fortran  program,  since  its 
last  line  is  a  call  to  FBACK.  For  more  information,  see  the 
DO  Fortran  IV  User's  Manual,  chapter  II-5. 
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C  PROGRAM  CMULT  -  DG  FORTRAN  3  -  CAPT  ROBERT  WELLS,  DEC  1983 

C  THIS  PROGRAM  PERFORMS  COMPLEX  CONJUGATE  MULTIPLICATION  OF  TWO 
C  236  X  236  COMPLEX  FILES  ON  A  POINT  BY  POINT  BASIS. 

INTEGER  1 1FLNM< 7) , I2FLNM(7).  0FLNM(7> 

INTEGER  MAIN(7) ,  MS (2) ,  SI ( 2 ) > S2 <  2 ) , S3 (2 ) 

COMPLEX  INI (2048) ,  IN2 < 2048 ), OUT (2048) 

CALL  I0F<3, MAIN.  I1FLNM,  I2FLNM, OFLNM. MS. SI. S2.  S3) 

CALL  0PEN(1.  1 1FLNM.  1,  IER  ) 

CALL  CHECK (IER) 

CALL  OPEN (2, I2FLNM.  1,  IER) 

CALL  CHECK (IER) 

CALL  OPEN (3. OFLNM, 3,  IER) 

CALL  CHECK (IER) 

DO  3  1-0.  31 

CALL  ROBLKd.  32*1,  INI,  32,  IER) 

CALL  CHECK (IER) 

CALL  RDBLK(2, 32*1,  IN2,  32,  IER) 

CALL  CHECK (IER) 

DO  2  J-l.  2048 

0UT(J)-IN1(J)*C0NJ0(IN2(J> > 

2  CONTINUE 

CALL  WRBLK(3.  32*1.  OUT, 32, IER) 

CALL  CHECK (IER) 

3  CONTINUE 

CALL  RESET 
STOP  CMULT 
END 
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C  PROGRAM  ADDCMP  -  DG  FORTRAN  -  CAPT  ROBERT  WELLS.  DEC  1983 

THIS  PROGRAM  ADDS  TWO  COMPLEX  FILES  TO  GET  ONE  COMPLEX  FILE. 

THE  FILES  MUST  BE  256  X  256  COMPLEX  FILES. 

INTEGER  AFLNM<  7) .  BFLNM<  7) .  CFLNM<  7) . MAIN! 7) , MS(2) > SI <2 ) . S2(2>  >  S3 (2) 
COMPLEX  AMAT <  256.  4 ) .  BMAT <  256.  4 ) .  CMAT ( 256. 4 ) 

CALL  I0F(3< MAIN.  AFLNM,  BFLNM.  CFLNM.  MS. SI.  S2,  S3) 

CALL  OPEN(l, AFLNM,  1,  IER ) 

CALL  CHECK! IER) 

CALL  OPEN (2.  BFLNM.  1. IER) 

CALL  CHECK! IER) 

CALL  OPEN! 3, CFLNM. 3. IER) 

CALL  CHECK! IER) 

DO  2  1-0.63 

CALL  RDBLK! 1,  16*1.  AMAT,  16.  IER) 

CALL  CHECK! IER) 

CALL  RDBLK! 2,  16*1.  BMAT,  16.  IER) 

CALL  CHECK! IER) 

DO  1  K-1.4 

DO  1  3-1.  256 

CMAT! J.  K ) -AMAT 1 J,  K)+BMAT! J,  K) 

1  CONTINUE 

CALL  URBLK!3.  16*1.  CMAT. 16.  IER) 

CALL  CHECK! IER) 

2  CONTINUE 
CALL  RESET 

STOP " <7><7><7><7> ADDCMP " 

END 
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C  PROGRAM  VDTOCP  -  DG  FORTRAN  5  -  CAPT  ROBERT  WELLS,  DEC  1983 

C  THIS  PROGRAM  CONVERTS  A  256  X  256  PACKED  VIDEO  FILE  TO  A 

C  256  X  256  COMPLEX  FILE 


i 

*  i 
_  # 

»  * 

4 

•1 


INTEGER  I FLNM ( 7 ) ,  OFLNM ( 7 > 

INTEGER  MAIN<7), F3<7)>  MS(2) , SI (2>, S2( 2) > S3 (2) 
INTEGER  TEMP < 512).  IN(2048> 

COMPLEX  OUT ( 2043 ) 


>5 

I 


CALL  I OF (2. MAIN,  I FLNM, OFLNM.  F3, MS.  SI, S2,  S3) 

CALL  OPENU,  I  FLNM,  1,  IER  ) 

CALL  CHECK (IER) 

CALL  OPEN <2, OFLNM, 3. IER) 

CALL  CHECK (IER) 

DO  3  1*0.  31 


CALL  RDBLK( 1,  2*1.  TEMP, 2, IER  > 
CALL  CHECK (IER) 

CALL  UNPACK<512,  TEMP,  IN) 


DO  2  J- 1,2048 

OUT( J)=CMPLX( FLOAT ( IN ( J) >.  0.  0> 
2  CONTINUE 


CALL  WRBLK(2. 32*1. OUT.  32. IER) 
CALL  CHECK (IER) 


3  CONTINUE 

CALL  RESET 
STOP  VDTOCP 
END 
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C  PROGRAM  CPTOVD  -  DG  FORTRAN  3  -  CAPT  ROBEKT  WELLS,  DEC  1983 

THIS  PROGRAM  FINDS  THE  MAXIMUM  AND  MINIMUM  OF  A  256  X  236 
COMPLEX  FILE  AND  GIVES  THE  LOCATION  OF  THE  MAXIMUM.  ALSO.  THE 
PROGRAM  CONVERTS  THE  FILE  TO  PACKED  VIDEO  FORMAT  BY  PERFORMING 
A  LINEAR  SCALING  TO  THE  RANGE  OF  0  -  15  AND  TRUNCATING. 

INTEGER  I FLNM < 7 ) .  OFLNM <  7 ) 

INTEGER  MAIN< 7) . F3 ( 7 ) . MS<2>, SI (2), S2<2).  S3<2) 

INTEGER  TEMP ( 512) >  OUT ( 2048 > ,  IMAX, JMAX, MAXCOL. MAXROW 
COMPLEX  IN <2048) 

REAL  RMAC. MAXVAL. MINVAL. RANGE 

CALL  I0F<2. MAIN.  IFLNM.  OFLNM.  F3.  MS,  SI.  S2.  S3) 

CALL  0PEN<1. IFLNM. 1. IER) 

CALL  CHECK (IER) 

CALL  0PEN<2, OFLNM, 3. IER) 

CALL  CHECK (IER) 

MAXVAL-O.  0 
MINVAL-1E30 

DO  3  1-0. 31  ;  FIND  MAX  AND  MIN 

CALL  RDBLK< 1, 32*1. IN. 32. IER) 

CALL  CHECK (IER) 

DO  3  J-  1.2048 

RMAC-CABS< IN< J) ) 

IF(RMAO.  LE.  MAXVAL >00  TO  2 
MAXVAL-RMAC 
IMAX**  I 
JMAX-J 

2  IF(RMAO.  LT.  MINVAL )MINVAL«RMAC 

3  CONTINUE 

RANGE-MAX VAL-M I NVAL 

MAXCOL-MOD <  JMAX- 1 . 236  >  + 1 
MAXR0W-8*IMAX+INT<  < JMAX-1 ) /236 ) *1 

TYPE-MAX  COL  -  ".MAXCOL  i  REPORT  MAX  IMFO. 

TYPE-MAX  ROW  -  ".MAXROW 

TYPE-MAX  VAL  -  ".MAXVAL 
TYPE-MIN  VAL  -  ".MINVAL 

DO  3  1-0,  31 

CALL  RDBLH< 1. 32*1. IN. 32. IER) 

CALL  CHECK (IER) 

DO  4  J- 1.2048 

OUT <  J)-ANINT (13.  0*(CADS< IN< J ) > -MINVAL > /RANGE ) 

4  CONTINUE 

CALL  REPACK <31 2. OUT, TEMP) 

CALL  WRBLK<2. 2*1. TEMP,  2.  IER) 

CALL  CHECK (IER) 

3  CONTINUE 

CALL  RESET 

STOP " <7><7><7><7>CPT0VD" 

END 
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PROGRAM  STATISTICS  -  DO  FORTRAN  9  -  CAPT  ROBERT  WELLS.  DEC  1983 

THIS  PROGRAM  COMPUTES  THE  MEAN  AND  THE  MEAN  DEVIATION  FOR  A 
NEIGHBORHOOD  AROUND  EACH  PIXEL  OF  THE  INPUT  FILE  < PACKED  VIDEO 
FORMAT).  THE  NEIGHBORHOOD  SIZE  AND  SHAPE  ARE  USER  SELECTABLE. 
THE  PROGRAM  USES  A  FAST  ALGORITHM  FOR  COMPUTING  THE  MEAN. 


INTEOER  0FLNM<7>. AFLNM < 7 > . DFLNM < 7 > 

INTEGER  TEMP ( 913) • OR 10 (296. 32) .  AVG(296>  8).  DEV <296.  8) 
INTEOER  MASKS Z. SHAPE.  XOFF. YOFF. PIXROW. BUFROW.  COL.  ROW 
INTEGER  HSUH.  FCOLSM. NRMSUH. CHANGE 
REAL  MSKCNT 


ACCEPT-ENTER  ORIGINAL  FILENAME  ->" 
READ! 11. l)OFLNM(l) 

FORMAT <S13> 

ACCEPT "ENTER  AVERAGE  FILENAME  ->" 
READ  <11. 1 ) AFLNM < 1 > 

ACCEPT "ENTER  DEVIATION  FILENAME 
R£AD< 11.  1 > DFLNM < 1 ) 


GET  FILENAMES 


2  ACCEPT "ENTER  MASK  SIZE  ->".MA8KSZ  :  GET  MASKSIZE  AND  SHAPE 
IF<  (MASKSZ.  NE.  3).  AND.  <MASKSZ.  NE.  9).  AND.  (MASKSZ.  NE.  7).  AND. 

•  < MASKSZ-  NE.  11).  AND.  < MASKSZ.  NE.  17).  AND.  (MASKSZ.  NE.  29))00  TO  2 

3  ACCEPT "CHOOSE  MASK  SHAPE:  O-HOR.  1-VERT.  2-SQUARE  ->". SHAPE 
IF <  (SHAPE.  NE.  0).  AND.  (SHAPE.  NE.  1).  AND.  (SHAPE.  NE.  2>)00  TO  3 


CALL  OPEN( 1.  OFLNM.  1.  IER  > 
CALL  CHECK < IER) 

CALL  OPEN (2. AFLNM.  3.  IER) 
CALL  CHECK (IER) 

CALL  OPEN < 3. DFLNM. 3.  IER) 
CALL  CHECK (IER) 


:  OPEN  FILES 


IF (SHAPE.  EQ.  0)00  TO  4 
IF  (SHAPE.  EQ.  1)00  TO  9 
XOFF-( MASKSZ- l>/2 
YOFF-XOFF 
00  TO  6 

4  XOFF* < MASKSZ— 1 )/2 
YOFF-O 
00  TO  6 
9  X OFF-0 

YOFF* ( MASKSZ -1 )/2 
6  CONTINUE 


)  FIND  X  ti  Y  OFFSETS 


M8KCNT*<2*X0FF+1)*(2*Y0FF+1)  <  NO.  OF  PIXELS  IN  MASK 

DO  7  1-1.912  <  FILL  IN  A  VO  FILE  BORDERS.  TOP  !■  BOTTOM 

TEMP < 1 ) -69939 
7  CONTINUE 

CALL  WRBLX(2<  0.  TEMP.  2.  IER) 

CALL  CHECK (IER) 

CALL  URBLK<2.  2.  TEMP.  1.  IER) 

CALL  CHECK! IER) 

CALL  WRBLK(2. 61.  TEMP. 2.  IER) 

CALL  CHECK (IER) 

CALL  WR1UM2.  63.  TEMP.  1.  IER) 

CALL  CHECK (IER) 


v> 


PACE 


2 


00  8  1-1.512  i  FILL.  IN  DEV  FILE  BORDERS.  TOP  It  BOTTOM 
TEMP  ( I  )— O 
8  CONTINUE 


CALL  WRBLKO.  0.  TEMP.  2.  IER ) 
CALL  CHECK (IER) 

CALL  WRBLHO.  2.  TEMP.  1.  IER) 
CALL  CHECK (IER) 

CALL  WRBLK ( 3. 61 >  TEMP . 2.  IER) 
CALL  CHECK (IER) 

CALL  WRBLKO.  A3.  TEMP,  1.  IER) 
CALL  CHECK (IER) 


DO  9  1-0.3  (READ  IN  FIRST  8  BLOCKS  OF  ORIC 
CALL  RDBLK( 1,2*1.  TEMP,  2,  IER) 

CALL  CHECK (IER) 

ROU— 8*1+1 

CALL  UNPACK<  912.  TEMP.  ORlOd,  ROW)) 

9  CONTINUE 


MSUM-0  i  COMPUTE  FIRST  AVERAGING  VARIABLES 
DO  10  J— 0, 2*Y0PP 
ROW— 1 3— YQFF+ J 
DO  10  K— O,  2*X0FF 
COL— 1 3— XOFF+K 
H8UM-MBUH+0R 10 ( COL,  ROW) 

10  CONTINUE 

CHANOE-O 
ROW— 1 3+ YOFF 
DO  11  K— 0, 2*X0FF 
COL— 13- XOFF+K 

CHANGE— CHANGE +OR 10 (COL. ROW) 

11  CONTINUE 


NRMSUM-M8UM-CHANGE 


DO  13  1-1.28  i  COUNT  OF  BUFFER  LOADS 

DO  12  BUFROW— 1,  8  i  ROW  COUNT  WITHIN  2K  BUFFERS 
P I XROW-8* I +BUFR0W+4 

CALL  DOROW(ORIG.  A VC,  DEV. PIXROW, BUFROW, XOFF,  YOFF.  MSKCNT. 
•  MSUM, NRMSUM) 

TYPE "ROW  NUMBER: ".PIXROW 


12  CONTINUE 

CALL  REPACK ( 912. AVO, TEMP)  i  WRITE  2K  BUFFERS  OF  RESULTS 
CALL  WRBLK(2. 2*1+1.  TEMP.  2. IER) 

CALL  CHECK (IER) 

CALL  REPACK (912.  DEV. TEMP) 

CALL  WRBLKO.  2*1+1.  TEMP,  2.  IER) 

CALL  CHECK! IER) 


CALL  RDBLKd.  2*I+«.  TEMP.  2,  IER)  »  READ  IN  NEXT  2K  OF  ORIC 
CALL  CHECK (IFR) 

ROW-MOD (8* I +24.  32)+l 
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CALL  UNPACK< 912.  TEMP.  0RI0<1.  ROW) ) 

13  CONTINUE 

00  14  8UFR0W»1.8  i  DO  LAST  2K  BUFFER 
P I XR0W-236+BUFR0W 

CALL  OOROUCORIO.  AVO. DEV. PIXROM.  BUFROU. XOFF.  YQFF. MSKCNT. 

•  MSUM.  NRMSUM) 

TYPE "ROW  NUMBER: ".PIXROU 

14  CONTINUE 

CALL  REPACK <312.  AVO. TEMP) 

CALL  WR8LK < 2.  39. TEMP.  2.  IER  > 

CALL  CHECK < IER) 

CALL  REPACK <312. DEV.  TEMP) 

CALL  UR8LKO.  39.  TEMP.  2.  IER) 

CALL  CHECK < IER) 

CALL  RESET 
STOP  STATS 
END 


SUBROUTINE  DOROWtORIG,  AVG,  DEV,  PIXROW. BUFROU. XOFF,  YOFF.  MSKCNT 
•  I1SUM.  NRMSUH) 

I NTEOER  TEMP  <  5 1 2  > .  OR  I Q  <  256.  32  ) ,  AVO  <  256.  8  > .  DEV  ( 256.  8  ) 

INTEGER  XOFF.  YOFF,  PIXROW.  BUFROW.  ROW.  COL 
INTEGER  HSUH, FCOLSM,  NRMSUH, CHANGE. PI XCOL. DSUH 
REAL  MSKCNT 

DO  20  K-l.  12  j  FILL  IN  FIRST  12  PIXELS 
AVG(K.  BUFROW) “15 
DEV(K> BUFROW ) -O 

20  CONTINUE 

DO  21  K-245.  256  »  FILL  IN  LAST  12  PIXELS 
AVGiK, BUFROW) -15 
DEV <K.  BUFROW >-0 

21  CONTINUE 

CHANGE— O 

ROW-MOD  (PIXROW-*- YOFF- 1.  32)  +  l 
DO  22  K— O. 2*X0FF 
COL— 13-XOFF+K 

CHANGE— CHANGE+OR 10 < COL.  ROW) 

22  CONTINUE 


MSUH-NRMSUH+CHANOE 


i  NEW  MATRIX  SUM 


CHANCE— 0 

ROW-HOD < P I XROW-YOFF-1 .  33>+l 
DO  23  K-O. 2* XOFF 
C0L-13-X0FF+K 

CHANGE— CHANOE+OR 1 0 ( COL,  ROW) 
23  CONTINUE 


NR HSUH— MS UM-C HA NOE 


i  NEW  NEXT.  ROW.  MSUH 


FCOLSM— O  i  COMPUTE  NEW  FIRST  COL  SUM 

COL- 13- XOFF 
DO  24  U-0.  2*Y0FF 

ROW-MOD ( P I XROW-YOFF+U- 1 .  32 ) ♦! 

FCOLSM— FCQLSM+OR 10 ( COL,  ROW) 

24  CONTINUE 


AVO( 13.  BUFROW )«ANINT(MSUM/M8KCNT)  i  PIXEL  13  AVERAGE 

DO  26  PIXCOL-14. 244  »  AVO  ROUTINE  FOR  PIXELS  14-244 

CHANGE-0 
COL-P I XCOL+XOFF 
DO  29  J— 0, 2*Y0FF 

ROW  MOD  <  P I XROW— YOFF»V- 1 ,  32)«-l 
CHANGE -CHANGE +0R I G < COL,  ROW) 

CONTINUE 

rmw  i"i  wuuan 

AVO ( P I XCOL. BUFROW )-ANINT(MSUM/MSKCNT) 

FCOLSM— 0 
COL-P I XCOL-XOFF 


PACE 


2 


DO  26  U-O, 2*Y0FF 

ROW“MOD <  P I XROW-YOFF  *  J-i » 32)+I 
FC0LSM»FC0L3M+0R IC ( COL.  ROW) 

26  CONTINUE 

00  28  PIXC0L*13,244  t  DEVIATION  ROUTINE 
DSUH*© 

00  27  J-O.2*Y0FF 

ROW-MOD<PIXROW-YOFF+J-I.  32)+l 
00  27  K-0.2#X0FF 

COL-P I XCOL-XOFF+K 

DSUN-DSUH+ABS (OR  10 (COL. ROW) -AVO (PIXCOL.  BUFROW ) ) 

27  CONTINUE 

DEV( PIXCOL, BUFROW) «AN I NT  <  4*D8UM/MSKCNT ) 
IF(DEV(PIXCOL. BUFROW)  CT.  19)DEV(PIXC0L.  BUFR0U)-19 

28  CONTINUE 

RETURN 

END 


C  PROGRAM  KEDGE  -  DO  FORTRAN  5  -  CAPT  ROBERT  WELLS.  DEC  1983 

C  THIS  PROGRAM  EDOES  THE  INPUT  FILE  (PACKED  VIDEO  FORMAT)  WITH  THE 

C  KIRSCH  OPERATOR.  THE  COEFFICIENTS  ARE  NOT  3  AND  -3  BUT  THE  RATIO 

C  IS  3  TO  -3.  THE  COEFFICIENTS  WERE  REDUCED  TO  KEEP  THE  EDGE  OUTPUT 

C  IN  THE  0-19  RANGE.  THE  OUTPUT  IS  396  X  396  PACKED  VIDEO  DATA. 

INTEGER  INMAT(396.  32).  OUTMAT (336. 16) . TEMP ( 1034). TOGGLE 
INTEGER  IFLNM<7). 0FLNM<7> 

INTEGER  MAIN(7>.  F3(7>.MS(2)»S1<2).S3<2).S3<2) 

CALL  10Pt3.MAlN.lFLNM.0FLNM.F3.M8.  SI.  S3.  S3) 

CALL  OFENtl,  IFLNM.  1.  IER) 

CALL  CHECK (IER) 

CALL  OPEN (3. OFLNM.  3.  IER) 

CALL  CHECK (IER) 

DO  1  1*1.912  i  FILL  FIRST  AND  LA8T  8  ROUS  OF 
TEMP < I >-0  j  OP  RESULTS  WITH  ZERO 

1  CONTINUE 

CALL  URBLKO.O.  TEMP.  3.  IER) 

CALL  CHECK (IER) 

CALL  WRBLK < 2.  63.  TEMP.  3.  IER) 

CALL  CHECK (IER) 

CALL  RDBLKd.O.  TEMP.  4.  IER)  i  READ  IN  FIRST  16  R0W8 
CALL  CHECK (IER) 

CALL  UNPACK! 1024. TEMP.  INMAT) 

T000LE*0  i  FLAG  TO  SHOW  BUFFER  WRAP-AROUND 

DO  9  1*1.  13  »  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  RDBLK( 1.4*1.  TEMP.  4.  IER)  > READ  NEXT  16  ROUS 
CALL  CHECK (IER) 

CALL  UNPACK!  1024.  TEMP.  INMATd.  17-1 6*T0GGLE )  > 

DO  3  K-i.  16  i  COUNTER  FOR  OUTPUT  BUFFER  ROUS 

DO  2  J-l.S  i  ZERO  FILL  FIRST/LAST  B  COLUMNS 
OUTHAT ( J.  K ) *0 
OUTMAT ( J+248.  K>*0 

2  CONTINUE 

DO  3  J-9.248  »  EDGE  THIS  ROW 

CALL  EDGER3< J.  K.  TOGGLE.  INHAT.  OUTMAT) 

3  CONTINUE 


CALL  REPACK (1034.  OUTMAT.  TEMP) 

CALL  WRBLK (3,  4*1-2. TEMP. 4, IER)  *  WRITE  16  ROUS  OF  RESULTS 
CALL  CHECK ( IER) 

IP (TOGGLE.  EG.  1)00  TO  4 

TOGGLE* 1  <  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 

GO  TO  9 

4  TOOOLE-O 

9  CONTINUE 
CALL  RESET 

STOP  -  <7><7><7X7>KED0E  " 

END 


SUBROUTINE  EDCER3< J. K, TOGGLE,  INMAT. QUTMAT) 

INTEGER  J. K. TOGGLE.  INMAT (256, 32) . 0UTMAT<236.  16) 

INTEGER  Jl. J2,  J3,  Kl.  K2,  K3 
REAL  A.  B.  C.  O.  E.  F.  0,  H.  MAXIM 

IF (TOGGLE. EG.  1)00  TO  1 
Kl-K+7 
K2-H+B 
K3-K+9 
00  TO  2 

1  Kl-M0D(K+22. 32>*1 
K2-M0D ( K+23,  32>>1 
K3-M0D(K*24.  32)*1 

2  Jl-J-1 
J2-J 
J3-J+1 

A-0.  667* ( INMAT ( Jl •  Kl ) ♦INMAT ( 02,  Kl  )♦ INMAT ( J3»  HI ) )  — 

*  0.  400*( INMAT (J3. K2)+TNMAT( J3,  K3)+INMAT( J2<  K3>* 

*  INMAT ( Jl ,  K3) ♦ INMAT ( Jl, K2) ) 

B-0  667* ( INMAT (J2,  Kl )*INMAT< J3.  Kl >+INMAT< J3.  K2) )- 

*  0.  400*( INMAT( J3,  K3)+INMAT( J2<  K3)+INMAT( Jl.  K3)* 

*  INMAT (Jl, K2) + INMAT (Jl. Kl ) > 

C-0  667* ( INMAT (  J3.  Kl  )*INMAT(  J3,  K2)  +  INMAT ( J3,  K3>  >- 

*  O.  400*<  INMAT<  J2.  K3 ) +INMAT ( Jl .  K3)*INMAT< Jl »  K2)  + 

*  INMAT (Jl,  Kl >+INHAT( J2»  Kl ) ) 

D-0.  667* ( INMAT( J3.  K2»*INMAT< J3,  K3>*INMAT( J2.  K3> )- 

*  0.  400* ( INMAT (Jl, K31+INMAT ( Jl,  K2)*INMAT ( Jl . Kl ) ♦ 

*  INMAT (J2,  Kl )*INMAT( J3.  Kl ) ) 

E-0.  667* ( INMAT ( J3, K3 ) ♦ INMAT ( J2, K3)+INMAT( Jl,  K3) ) - 

*  O.  400* ( INMAT (Jl,  K2)+INMAT( Jl, Kl )+INMAT( J2,  Kl  )♦ 

*  INMAT ( J3,  Kl )♦ INMAT (J3. K2) ) 

F-0.  667* ( INMAT( J2.  K3)*INMAT( Jl.  K3)+INMAT( Jl. K2) >- 

*  0.  400*( INMAT( Jl,  Kl >*INMAT( J2,  Kl )+INMAT( J3,  Kl >♦ 

*  INMAT ( J3.  K2>* INMAT (J3.K3)) 

0-0.  667* ( INMAT (Jl, K3)+INHAT( Jl, K2)+INMAT( Jl, Kl>  >— 

*  0.  400*( INMAT( J2.  Kl )*INMAT( J3#  Kl )+INMAT( J3,  K2)* 

*  INMAT (J3, K3)+INMAT( J2.  K3>  > 

H— O.  667* ( INMAT( Jl,  K2)*INMAT( Jl,  Kl )+INMAT( J2,  Kl) )- 

*  0.  400*< INMAT (J3. Kl )*INHAT( J3.  K2)*INMAT< J3.  K3>* 

*  1NHAT( J2,  K3)*INMAT( Jl.  K3> ) 

MAXIM-AMAX1 ( A. B,  C.  0.  E. F.  0.  H> 

OOTMAT( J.  K)-ANINT( MAXIM) 

IF(OUTMAT(J. K).  OT.  13)0UTMAT(J.  K)«19  j  CHECK  FOR  CLIPPING 


RETURN 

ENO 


nonnn 


PACE 


1 


C  PROGRAM  WEDGE  -  DG  FORTRAN  9  -  CAPT  ROBERT  WELLS.  DEC  1963 

THIS  PROORAH  EDGES  THE  INPUT  FILE  (PACKED  VIDEO  FORMAT)  WITH 
A  MASK  OPERATOR  CALLED  THE  WEDGE  OPERATOR.  THE  OUTPUT  IS  SEPARATED 
INTO  FOUR  294  X  294  PACKED  VIDEO  FILES  ACCORDING  TO  WHICH  MASK 
PRODUCED  THE  MAXIMUM  VALUE.  IF  TWO  MASKS  OF  DIFFERENT  ORIENTATIONS 
PRODUCE  THE  MAXIMUM.  THAT  VALUE  GOES  TO  BOTH  OUTPUT  FILES. 

INTEGER  INMAT(294.  14). TEMP(912>. TOGGLE 

INTEOER  OMAT1 (294.  8) . 0MAT2(294. 8 ) .  0MAT3 ( 294.  8).  0MAT4(294.  8) 

INTEGER  INFLNM(7).0FLNM1(7).0FLNN2(7).QFLNM3(7).0FLNM4(7> 

ACCEPT "ENTER  INPUT  FILENAME 
READUl.  DINFLNMU) 

1  F0RMAT(S13) 

ACCEPT*ENTER  OUTPUT  FILENAME  #1 
READ (11.  l)OFLNHl(l) 

ACCEPT "ENTER  OUTPUT  FILENAME  *2 
READ (11.  1)0FLNM2(1> 

ACCEPT "ENTER  OUTPUT  FILENAME  43  ->" 

READ (11.  1 )0FLNH3< 1 ) 

ACCEPT "ENTER  OUTPUT  FILENAME  44  ->" 

READ (11.  1 )0FLNH4( 1 ) 

CALL  OPENO.  INFLNH.  I.  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 2.  OFLNM1.  3.  IER) 

CALL  CHECK (IER) 

CALL  OPENO.  0FLNM2.  3.  IER) 

CALL  CHECK! IER) 

CALL  OPEN (4.  0FLNM3.  3.  IER) 

CALL  CHECK! IER) 

CALL  OPENO.  0FLNH4.  3.  IER) 

CALL  CHECK! IER) 

DO  2  1-1.294  I  FILL  FIR8T  AND  LAST  4  ROWS  OF 
TEMP! I )— 0  »  OP  OUTPUT  FILES  WITH  ZERO 

2  CONTINUE 

CALL  WR8LK (2.  0.  TEMP.  1.  IER) 

CALL  CHECK! IER) 

CALL  WRBLK(2.  43.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLKO.  0.  TEMP.  1.  IER) 

CALL  CHECK! IER) 

CALL  WRBLKO.  43.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLKO.  0.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLKO.  43.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLKO.  0.  TEMP.  1.  IER) 

CALL  CHECK! IER) 

CALL  WRBLKO.  43.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  RDSLK( 1. 0. TEMP. 2.  IER)  > READ  IN  FIRST  8  ROWS 
CALL  CHECK (IER) 

CALL  UNPACK! 912.  TEMP.  INMAT) 

TOOCLE-O  i  FLAG  TO  SHOW  BUFFER  WRAP-AROUND 
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00  9  1*1.31  I  COUNTER  FOR  2M  BUFFER  LOADS 

CALL  RDBLK( 1.2*1.  TEMP. 2. IER>  i  READ  NEXT  S  ROWS 
CALL  CHECK (XER) 

CALL  UNPACK012.  TEMP.  INMAT! 1. 9-8*T000LE) > 

DO  3  K-1.8  i  COUNTER  FOR  OUTPUT  BUFFER  R0W8 

CALL  WR0W1 <K.  TOC OLE, INMAT.  0HAT1.  0MAT2.  0MAT3.  0MAT4) 

3  CONTINUE 

CALL  REPACK! 912. 0MAT1. TEMP) 

CALL  WRBLK(2. 2*1-1. TEMP. 2. IER)  »  WRITE  S  ROWS  OP  RESULTS 
CALL  CHECK! XER) 

CALL  REPACK ! 912.  0MAT2.  TEMP ) 

CALL  WRBLK!3. 2*1-1. TEMP.  2.  IER) 

CALL  CHECK! XER) 

CALL  REPACK ! 91 2.  0MAT3.  TEMP) 

CALL  WRBLK!4.  2*1-1.  TEMP.  2.  IER) 

CALL  CHECK! XER) 

CALL  REPACK ! 912. 0MAT4.  TEMP) 

CALL  WRBLK!9. 2*1-1. TEMP. 2.  IER ) 

CALL  CHECK! IER) 

IFiTOQOLE.  EQ.  1)00  TO  4 

TOOOLE-1  i  TOOOLE  THE  BUFFER  WRAP-AROUND  FLAO 
00  TO  9 

4  TOQOLE-O 
9  CONTINUE 

CALL  RESET 

STOP  "  <7><7X7X7>WEDC£  ■ 

END 


PACE 
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SUBROUTINE  WROU1 (K.  TOOOLE.  INMAT.  OMAT1,  0MAT2.  0MAT3.  0MAT4) 
INTEOER  K. TOOOLE.  INMAT (296. 16) . HEDGE 

INTEOER  OMAT 1(296.  8) .  QMAT2<296.  8 > .  0MAT3 ( 296.  8) . 0MAT4(296. 8) 
REAL  Al.  A2.  A3.  A4.  Bl.  B2.  B3.  B4.  Cl.  C2.  C3.  C4,  01.  02.  03.  04.  MAXIM 

IF ( TOOOLE.  EQ.  1)00  TO  1  j  FIND  INPUT  BUFFER  INDEXES 
K1-K4-3 
K2-H+4 
K3-K4-5 
00  TO  2 

1  K1-MOO(K*10.  16)4-1 
K2-MOO(K-*-ll.  16)4-1 
K3«M0D(K+12.  16)4-1 

2  CONTINUE 

DO  3  J-1.4  i  ZERO  FILL  FIRST/LAST  4  COLUMNS  OF  RESULTS 
OMATKd.  K)-0 
OMAT  1  ( J4-232.  K )  «0 
0MAT2(J.K)-O 
0MAT2(  J4-292.  K)-0 
0MAT3(J.K)-0 
0MAT3  <  J+292.  K  >  -0 
0MAT4(J.K)-0 
0MAT4(  J4-292.  K)-0 

3  CONTINUE 

00  4  J»9. 292  i  EDGE  THESE  COLUMNS 

J1«J~1 

J2-J 

J3-J4-1 

Al-O.  600*  < INMAT(  Jl.  Kl  )4-INMAT(  J2,  K1  )4-lNHAT(  J3.  KI  )  >- 
4-  O.  300»(INMAT(dl,K2)4-INMAT(v)2.  K2)4-INMAT(J3.  K2)4- 

4-  INMAT(  Jl.  K3)4-INMAT(  J2,  K3)4-INMAT(  J3.  K3) ) 

A2-0.  600*(  INMAT(J1.  K3)4-INMAT(  J2.  K3)4-INMAT(  J3.  K3)  >- 
4-  0.  300*(  INMAT  (Jl.  K2)4-INMAT(  J2,  K2)4-INMAT(  J3.  K2)4- 

♦  INMAT(  Jl.  K1)4-INMAT(  J2.  Hi  )4-INMAT(  J3.  Kl  )  ) 

A3— Al 
A4— A2 

Bl-O.  600* ( INMAT ( J2.  Kl  >-*-INMAT(  J3.  Kl  )4-INMAT(  J3.  K2)  )- 
4-  0.  300*(INMAT(J1.K1)4-INMAT(J2.  K2)4-INMAT(  J3.  K3>4- 

♦  INMAT <  Jl.  K2)4-INMAT(  Jl.  K3 )  + INMAT (J2>  K3)  ) 

B2-0.  600* (  INHAT<  Jl.  K2)4-INMAT(  Jl.  K3)4-INMAT(  J2.  K3)  >- 

♦  0.  300*  ( INMAT  (Jl.  Kl  )+INMAT( J2<  K2  )♦  INMAT  (J3<  K3)4> 

♦  INMAT (J2.  Kl  )4>INMAT(  J3.  Kl  >*INMAT(  J3.  K2)  ) 

B3— II 

i4»-n 

Cl-O.  600* ( INMAT (J3.  Kl  )4-INMAT(  J3,  K2)4-INMAT(  J3.  K3>  )- 
4-  0.  300*(INMAT(J2.Ki>4-INMAT(J2.K2)4-INMAT(J2.K3)4. 

♦  INMAT(  Jl.  Kl  )♦ INMAT (Jl.  K2 )+ INMAT (Jl«  K3)  > 


♦ 

♦ 


C2-0.  600* < INMAT ( Jl>  K1 )+INMAT< Jl.  K2 )♦ INMAT <  Jl » K3) )- 
0.  300* ( INMAT ( J2>  HI )  +  INMAT  ( J2,  K2)  INMAT  <  J2,  K3>* 

INMAT (J3.  K1 )♦ INMAT ( J3» K2 ) ♦ I NMAT ( J3. K3) > 

C3— Cl 
C4— C2 

01-0.  600* < INHAT< Jl.  K2 )  +  INMAT ( Jl .  Kl >*INMAT< J2.  Kl) >- 
0.  300* <  INMAT ( Jl .  K3)*INMAT ( J2<  K2)*INMAT<  J3.  Kl  >•*■ 

*  INMAT<J2<  K3 ) +INMAT ( J3< K3)*INMAT< J3< K2) > 

02*0.  600* < INMAT <J2< K3)+1NMAT  <  J3< K3)*INMAT< J3.  K2) >- 

*  0.  300* ( INMAT< Jl .  K3 ) ♦ INMAT <  J2<  H2 ) +INMAT <  J3«  Kl ) ♦ 

*  INFIAT ( Jl <  K2 ) +INMAT <  Jl <  Kl  )«■  INF<AT<J2<  Kl)  ) 

D3-- D1 
04—02 

MAXIM-AMAX1 <A1.  A2< A3.  A4>  Bl.  B2<  B3.  B4.  Cl.  C2.  C3.  C4,  01.  D2<  03. D4) 

0MAT1 < J.  K)»0 
0MAT2<J.  K)-0 
0MAT3< J.  K)«0 
0MAT4<J.  K)-0 

MEDOE-ANINT<MAXIM)  j  INTEGRAL  EDGE  VALUE 

IF<MED0E.  OT.  1S)MED0E-15  i  TEST  FOR  OVERFLOW 
IF<MEDCE. LT. 0)MEDOE-0  i  AND  UNDERFLOW 

IF<A1.  EQ.  MAXIM.  OR.  A2.  EQ.  MAXIM.  OR.  A3.  EQ.  MAXIM.  OR.  A4.  EQ.  MAXIM) 

*  0MAT1 < J. K) “HEDGE  i  ORIENTATION  41 

IF<B1.  EO.  MAXIM.  OR.  B2.  EO.  MAXIM.  OR.  B3.  EO.  MAXIM.  OR.  B4.  EO.  MAXIM) 

*  QMAT2<J. K) -HEDGE  J  ORIENTATION  42 

IF <C1.  EO.  MAXIM.  OR.  C2.  EO.  MAXIM.  OR.  C3.  EO.  MAXIM.  OR.  C4.  EO.  MAXIM) 

*  0MAT3< J. K)*MEOOE  j  ORIENTATION  43 

IF<01.  EO.  MAXIM.  OR.  02.  EO.  MAXIM.  OR.  03.  EQ.  MAXIM.  OR.  04.  EO.  MAXIM) 

*  0MAT4< J. K) -HEDGE  »  ORIENTATION  44 

4  CONTINUE 


RETURN 

END 


nonnnn 
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C  PROGRAM  TEDCE  -  00  FORTRAN  9  -  CAPT  ROBERT  WELLS.  DEC  1983 

THIS  PROGRAM  DETECTS  LINES  IN  THE  INPUT  FILE  (PACKED  VIDEO  FORMAT) 
WITH  A  MASK  OPERATOR  FOR  LINE  DETECTION  AT  FOUR  ORIENTATIONS. 

THE  OUTPUT  IS  SEPARATED  INTO  FOUR  296  X  296  PACKED  VIDEO  FORMAT 
FILES  ACCORDING  TO  WHICH  MASK  PRODUCED  THE  MAXIMUM  VALUE.  IF  TWO  OR 
MORE  MASKS  PRODUCE  THE  SAME  MAXIMUM.  THE  OUTPUT  GOES  TO  EACH 
RESPECTIVE  FILE. 

INTEGER  INMAT(296.  16).  TEMP(912><  TOGGLE 

INTEGER  OMAT 1 ( 296.  8 ) .  0MAT2 ( 296.  8 ) . 0MAT3 ( 296.  8  > .  0MAT4 ( 296. 8  > 

INTEGER  INFLNM<7).  0FLNM1 (7).  0FLNM2<7>. 0FLNM3<7>.  0FLNM4(7> 

ACCEPT "ENTER  INPUT  FILENAME 
READ! It.  1 > INFLNH! 1 ) 

1  FORMAT (SI 3) 

ACCEPT-ENTER  OUTPUT  FILENAME  41 
READ (11. l)OFLNMl(l) 

ACCEPT-ENTER  OUTPUT  FILENAME  42  ->“ 

READ (11.  1>0FLNM2(1> 

ACCEPT-ENTER  OUTPUT  FILENAME  43  ->“ 

READ (11. 1 )0FLNM3( 1 ) 

ACCEPT-ENTER  OUTPUT  FILENAME  44  ->" 

READdl.  1  )0FLNH4(  1 ) 


CALL  CIPENO.  INFLNM.  1.  IER) 
CALL  CHECK! IER) 

CALL  0PEN(2. 0FLNM1.  3.  IER) 
CALL  CHECK (IER) 

CALL  OPEN! 3. 0FLNM2.  3.  IER) 
CALL  CHECK! IER) 

CALL  OPEN! 4. 0FLNM3.  3.  IER) 
CALL  CHECK! IER) 

CALL  OPEN ! 9.  0FLNM4.  3.  IER) 
CALL  CHECK (IER) 


DO  2  1-1.296  i  FILL  FIRST  AND  LAST  4  ROUS  OF 
TEMP ( I )— O  i  OP  OUTPUT  FILES  WITH  ZERO 

2  CONTINUE 

CALL  UR8LK!2. 0. TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLK ( 2.  63.  TEMP.  I.  IER ) 

CALL  CHECK! IER) 

CALL  URBLKO.O.  TEMP,  t.  IER) 

CALL  CHECK (IER) 

CALL  WRBLK (3, 63.  TEMP.  1,  IER) 

CALL  CHECK ( IER) 

CALL  WRBLK(4.0.  TEMP,  I.  IER) 

CALL  CHECK! IER) 

CALL  WRBLK (4. 63. TEMP.  1,  IER) 

CALL  CHECK (IER) 

CALL  WRBLKO.O.  TEMP.  1.  IER) 

CALL  CHECK (IER) 

CALL  WRBLK ( 9.  63. TEMP.  1,  IER) 

CALL  CHECK! IER) 


CALL  RDBLKd.O.  TEMP.  2.  IER)  i  READ  IN  FIRST  8  ROWS 
CALL  CHECK (IER) 

CALL  UNPACK(912,  TEMP,  INMAT) 
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TOCCLE-O  i  FLAG  TO  SHOW  SUFFER  WRAP-AROUND 

DO  S  1*1.31  t  COUNTER  FOR  2K  BUFFER  LOADS 

CALL  RDBLK< 1.2*1.  TEMP, 2.  IER>  >  READ  NEXT  8  ROWS 
CALL  CHECK < IER ) 

CALL  UNPACK < 912.  TEMP. INHAT< 1 . 9-8*T0CCLE> ) 

DO  3  K"l. 9  i  COUNTER  FOR  OUTPUT  BUFFER  ROUS 

CALL  TROW1 (K. TOGGLE.  INHAT, OMAT1. 0MAT2.  0HAT3,  0HAT4) 

3  CONTINUE 

CALL  REPACK (912.  OMAT1 . TEMP > 

CALL  WRBLK<2.  2*1-1. TEMP. 2, IER)  ;  WRITE  8  ROUS  OF  RESULTS 
CALL  CHECK! IER) 

CALL  REPACK (912. 0HAT2,  TEMP) 

CALL  WR8LKO,  2*1-1,  TEMP.  2,  IER) 

CALL  CHECK (IER) 

CALL  REPACK ( 912. 0MAT3. TEMP) 

CALL  URBLK(4,  2*1-1,  TEMP. 2.  IER) 

CALL  CHECK! IER) 

CALL  REPACK( 912. 0MAT4. TEMP) 

CALL  URBLKO.  2*1-1.  TEMP.  2.  IER) 

CALL  CHECK! IER) 

IF(TOOOLE.  EQ.  1)00  TO  4 

TOGGLE* 1  J  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 
GO  TO  9 

4  TOOOLE-O 
9  CONTINUE 

CALL  RESET 

STOP  “  <7><7X7X7>TEDQE“ 


SUBROUTINE  TROW1 (K> TOGGLE/  INMAT.  OMATl,  0MAT2,  0MAT3. 0MAT4> 


INTEGER  K. TOGGLE.  INMAT(2S6.  16).  A,  B.  C. D.  MAXIM 

INTEGER  OMATl (296.  8) .  0MAT2(296.  8) .  0MAT3C256.  8 ) . 0MAT4 ( 296. 8) 

IF (TOGGLE.  EG.  1) GO  TO  1  »  FIND  INPUT  BUFFER  INDEXES 

Kl-K+3 
K2-H+4 
K3-K+9 
GO  TO  2 

1  K1«M0D(K+10.  16)+1 
K2-M0D(K+ll.  16) >1 
K3*M0D ( K+ 12.  16)  +1 

2  CONTINUE 

DO  3  J»1.4  i  ZERO  FILL  FIRST/LAST  4  COLUMNS  OF  RESULTS 
OMATKJ.  K)-0 
OMATl (J+292.  K)-0 
0MAT2(J.  K)-0 
0MAT2( J+292.  K ) 

0MAT3(  J.  K>**0 
0MAT3( J+292.  K)-0 
0MAT4(J.  K)-0 
0MAT4 ( J+292.  H>«0 

3  CONTINUE 

DO  4  J-9.  292  i  SEP ER ATE  THE  TEMPLATE  EDGES 

Jl-J-1 

J2-J 

J3-J+1 

OMATl  (J.K)-O 
QHAT2 ( J.  K  )  *0 
0MAT3(J.K>-0 
0MAT4(J.  K)«0 

IF( INMAT ( J2.  K2>  .  EQ.  0)00  TO  4  :  CENTER  PIXEL  SHOULD  NOT  EQUAL  0 
A-INMAT( Jl.  K2) +INMAT ( J2.  K2)+INMAT( J3. K2) 

B-INMAT( Jl.  K1 )+INMAT( J2.  K2)+INMAT( J3. K3) 

C« INMAT ( J2.  K1 )  +  INMAT ( J2.  K2) +INMAT( J2. K3) 

D- INMAT ( Jl.  K3)+INMAT( J2.  K2)+INMAT( J3. K1 ) 

MAXIM-MAXO(A.B.C.D) 

IF  (MAXIM.  LT.  20)00  TO  4  'i  TEST  TO  ELIMINATE  NOISE 

IF(A.  EQ.  MAXIH>0MAT1<J.*>-19  *  ASS I ON  EDGE  VALUE  TO 

»  PROPER  OUTPUT  BUFFER (S) 

1F(B. EQ.  MAXIM)0MAT2( J.  H)-19 
IF(C.  EQ.  MAXIM)0MAT3( J. K)-19 
IF(D.  EQ.  MAXIM)0MAT4<  J. K)*19 

4  CONTINUE 

RETURN 

EM) 
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PROCRAM  HAM1  (SPATIAL  CLUSTER  RECOGNITION  ALGORITHM) 
(PARTIALLY  IMPLEMENTED  HERE) 

INTEGER  MAIN(7). IFLNM(7>. F2(7>, F3(7>. MS(2).  SI (2).  S2(2>.  S3 (2) 
REAL  XNMAT ( 128< 32) <  OUTMAT ( 128,  16) . ENPCT,  EDPCT.  ENTHR.  EDTHR 

CALL  IOF( 1. MAIN. IFLNM. F2. F3.  MS.  SI. S2.  S3) 

CALL  0PEN(1.  IFLNM.  1.  IER> 

CALL  CHECK (IER) 

CALL  OPEN (2,  "ENHC.  IR",  2.  IER) 

CALL  CHECK (IER) 

CALL  OPEN (3, “NTHR.  IR". 3, IER) 

CALL  CHECK (IER) 

CALL  OPEN (4.  "EDGE.  !R".2. IER) 

CALL  CHECK (IER) 

CALL  OPENO,  "DTHR.  IR".  3.  IER) 

CALL  CHECK (IER) 

CALL  OPEN ( 6. “CTHR.  IR". 2.  IER) 

CALL  CHECK (IER) 

CALL  0PEN(7. “CONN.  IR",  2. IER) 

CALL  CHECK (IER) 

ENHANCED  IMAGE  IS  FORMED  FROM  ORIGINAL  IMAGE 
CALL  NHANCEU.  2.  INMAT.  OUTMAT) 

ENPCT- 1.  0 

CALL  THRESH (2.  3. ENPCT.  ENTHR,  INMAT) 

EDOED  IMAGE  IS  FORMED  FROM  ENHANCED  IMAGE 
CALL  EG I MAO (2,  4. INMAT.  OUTMAT) 

EDPCT— 1 .  0 

CALL  THRESH (4,  S. EDPCT,  EDTHR. INMAT) 

EDOED  IMAGE  IS  REDUCED  BY  CONJUNCTIVE  THRESHOLDING 
CALL  C0NUUN(2. 4. 6. ENTHR. EDTHR.  INMAT, OUTMAT) 

EDOED  IMAGE  IS  REDUCED  BY  CONNECTEDNESS  TEST 
CALL  C0NNEC(6, 7. INMAT.  OUTMAT) 

CALL  RESET 

STOP " <7><7><7><7>HAM 1 " 


SUBROUTINE  NHANCE! INFIL, OUTFIL.  INMAT,  OUTMAT) 


r  ••• 


INTEGER  INFIL. OUTFIL,  TOGGLE 
REAL  INMAT< 128. 32).  OUTMAT <128,  16) 

DO  1  K-l.B  j  FILL  FIRST  AND  LAST  8  ROUS  OF 
DO  1  J-l.  128  i  RESULTS  WITH  ZERO 
OUTMAT (J,  K)«0.  0 

1  CONTINUE 

CALL  WRBLK (OUTFIL.  0.  OUTMAT. 8.  IER) 

CALL  CHECK! IER) 

CALL  WRBUUOUTFIL.  120.  OUTMAT.  8.  IER) 

CALL  CHECK! IER) 

CALL  RD9LK! INFIL.  0.  INMAT, 16.  IER)  i  READ  IN  FIRST  16  ROWS 
CALL  CHECK! IER) 

TOOOLE-O  i  FLAG  TO  SHOW  BUFFER  WRAP-AROUND 

DO  9  1-1.7  i  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  ROBLK! INFIL.  16*1.  INMAT! 1. 17-16»T000LE>.  16.  IER) 

CALL  CHECK! IER)  i  READ  IN  NEXT  16  ROWS 

DO  3  K— 1.  16  i  COUNTER  FOR  OUTPUT  BUFFER  ROWS 

DO  2  J-l.  8  i  ZERO  FILL  FIRST/LAST  8  COLUMNS 
OUTMAT !J.K)«0 
OUTMAT! Jf 120.  K)*0 

2  CONTINUE 

DO  3  J— 9,  120  J  ENHANCE  THIS  ROW 

CALL  NHANCl ! J. K. TOOCLE. INMAT. OUTMAT) 

3  CONTINUE 

CALL  WRBLK ! OUTFIL.  16*1 -B.  OUTMAT, 16. IER)  »  WRITE  16  ROWS 
CALL  CHECK! IER)  t  OF  RESULTS 

IF!TOQQL£.  EQ.  1)00  TO  4 

TOGGLE- 1  i  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 
00  TO  9 

4  TOOOLE-O 
9  CONTINUE 


TYPE "ENHANCED  IMAGE  COMPLETE* 

RETURN 

END 


SUBROUTINE  NHANC 1  ( J.  K, TOGGLE*  INMAT, OUTMAT ) 

INTEGER  J. K. TOGGLE*  Ji* J2*  J3. Kl. K2,  K3 

REAL  INMAT ( 128. 32).  OUTMAT <128.  16) , A. B. C, D. E. F.  O* H 

IF (TOGGLE.  EQ.  1)00  TO  1 
Kl-K+7 
K2-H+8 
K3-K+9 
00  TO  2 

1  Kl-M0D<K-t-22.  32)+l 
K2-M0D  ( K+23.  32 )  4*1 
K3-M0D(K«-24.  32>*1 

2  Jl-J-1 
J2-J 
J3-J+1 

A*10*< INMAT < Jl . Kl ) «■ INMAT ( J2, Kl ) +INMAT (J3. Kl )) - 

♦  ( INMAT ( J3. K2 )  + 1 NMAT ( J3.  K3 ) +INMAT( J2.  K3>* 

♦  INMAT ( Jl*  K3) ♦INMAT ( Jl.  K2 ) ) 

B»10*( INMAT (J2. Kl )+INMAT< J3. Kl >+INMAT( J3, K2> )- 
+  ( INMAT( J3. K3)  +  INMAT( J2.  K3 > +INMAT ( Jl .  K3)  + 

INMAT< Jl .  K2) +INMAT( Jl,  Kl > ) 

C-10*( INMAT (J3.  Kl ) ♦INMAT ( J3. K2)«-INMAT(  J3,  K3) >- 

♦  <  INMAT(J2.  K3)«-INMAT<J1.K3>  +INMAT (J1,K2)  + 

♦  INMAT ( Jl .  Kl )♦ INMAT (J2, Kl ) ) 

D-10* < INMAT < J3.  K2 ) ♦INMAT ( J3.  K3)«-INMAT<  J2,  K3> )- 

♦  <  INMAT  (Jl.  K3)+INMAT( Jl,  K2)+JNMAT( Jl,  Kl  )♦ 

♦  INMAT ( J2. Kl ) ♦INMAT ( J3,  Kl ) ) 

E-10*<INMAT(J3. K3)  + INMAT ( J2.  K3)«-INMAT<  Jl,  K3) )- 

♦  ( INMAT ( Jl. K2) +INMAT ( Jl, Kl )+INMAT <  J2,  Kl ) ♦ 

♦  INMAT (J3>  Kl )>INMAT ( J3. K2) ) 

F«10*( INMAT (J2.  K3)-*-lNMAT(Jl,K3)*INHAT<Jl,  K2)  )- 

♦  <  INMAT(  Jl.  Kl  )*1NMAT<  J2,  Kl  H-INMAT(  J3.  Kl  )♦ 

♦  INMAT (J3»  K2)>INMAT ( J3.  K3) ) 

0-l0*< INMAT( Jl.  K3 ) ♦! NMAT <J1.  K2)*INMAT< Jl, Kl ) )- 

♦  (INMAT<J2.  Kl)-*  INMAT  <J3.  K  l  >  •*■  I  NMAT  <  J3.  K2>«- 

♦  INMAT( J3. K3) ♦INMAT ( J2,  K3) ) 

H"10*< INMAT ( Jl , K2) +INMAT ( Jl . Kl ) +INMAT ( J2.  Kl ) )- 

♦  (INMAT( J3. Kl)fINMAT< J3.  K2 ) +INMAT( J3,  K3)  + 

+  INMAT( J2# K3)  +  INMAT<  Jl,  K3) > 

OUTMAT < J.  K)-AMAX1 (A.  B.  C,  D.  E,  F.  0.  H) 


D-l* 


RETURN 

END 


SUBROUTINE  THRESH< INFIL.  OUTFIL.  PCT.  THRLD.  INMAT) 
INTEGER  INFIL.  OUTFIL 

REAL  PCT, THRLD,  INMAT ( 128. 32).  RSUM, A VO, DEV 

RSUM-O.  0  i  GET  SET  TO  FIND  AVERAGE 

CALL  RDBLK! INFIL. 16. INHAT. 32. IER) 

CALL  CHECK! IER ) 

DO  1  K-3, 32  i  COUNTER  FOR  INPUT  BUFFER  ROWS 
DO  1  J-17.  112  »  COLUMN  COUNT 
RSUM-RSUH* 1 NHAT < J.  K ) 

1  CONTINUE 

CALL  RD8LK! INFIL.  48.  INHAT,  32.  IER) 

CALL  CHECK! IER) 

DO  2  K*l.  32  »  COUNTER  FOR  INPUT  BUFFER  ROWS 
DO  2  J-17, 112  t  COLUMN  COUNT 
RSUM-RSUH* 1 NHAT ! J.  K ) 

2  CONTINUE 

CALL  RDBLK! INFIL.  80.  INMAT.  32.  IER) 

CALL  CHECK! IER) 

DO  3  KM.  30  l  COUNTER  FOR  INPUT  BUFFER  ROWS 
DO  3  J-17. 112  i  COLUMN  COUNT 
RSUM-RSUH* 1 NHAT ! J.  K ) 

3  CONTINUE 

AVO-RSUM/ ! 96*92  > 

RSUM-O.  O  >  GET  SET  TO  FIND  DEVIATION 

CALL  RDBLK! INFIL. 16. INMAT. 32. IER) 

CALL  CHECK! IER) 

DO  4  K-3.  32  i  COUNTER  FOR  INPUT  BUFFER  ROWS 
DO  4  J-17. 112  i  COLUMN  COUNT 

R8UM-R8UH+  ( INMAT !  J.  K ) -AVQ ) **2 

4  CONTINUE 

CALL  RDBLK! INFIL.  48.  INHAT.  32.  IER) 

CALL  CHECK! IER) 

DO  3  K-1.32  »  COUNTER  FOR  INPUT  BUFFER  ROWS 

DO  3  J-17.  112  i  COLUMN  COUNT 

R8UM-R8UH+ ! INMAT ! J.  K l-AVO )**2 
3  CONTINUE 

CALL  RDBLK! INFIL.  80.  INHAT, 32,  IER) 

CALL  CHECK! IER) 

DO  6  K-1.30  i  COUNTER  FOR  INPUT  BUFFER  ROWS 
DO  6  J-17.  112  i  COLUMN  COUNT 

RSUM— RSUM* ! I NHAT I J,  K ) -AVG  >  **2 

6  CONTINUE 

DEV-SORT 1 RSUH/ < 96*92 ) ) 

THRLD— AVG *PCT*DEV 
TYPE "THRESHOLD  -  ",  THRLD 
00  8  1-0.3 

CALL  RDSUMINFIL.  32*1.  INMAT.  32,  IER) 

CALL  CHECK  HER) 

DO  7  K-1.32 

DO  7  J-l.  128 

IF!  INMAT!  J.  K).  LE.  THRLD)  INMAT!  J.  K)-0.  0 

7  CONTINUE 

CALL  WRSLK! OUTFIL,  32*1.  INHAT.  32.  IER) 

CALL  CHECK! IER) 

8  CONTINUE 


RETURN 

END 


SUBROUTINE  ESIMAC<INFIL<OUTFIL<  INMAT.  OUTMAT ) 


INTEGER  INFIL.  OUTFIL.  TOGGLE 
REAL  INMATU2B.  32).  OUTMAT (128.  16) 

00  1  K«l.  8  i  FILL  FIRST  AND  LAST  8  ROUS  OF 
DO  1  JM.  128  >  RESULTS  WITH  ZERO 

OUTMAT (J, K)-0.  0 

1  CONTINUE 

CALL  WRBLK (OUTFIL. 0. OUTMAT. 8. IER) 

CALL  CHECK! IER > 

CALL  WRBLK ! OUTF 1 L.  120. OUTMAT. 8. IER) 

CALL  CHECK! IER ) 

CALL  RDBLK! INFIL.  0.  INMAT.  16.  IER)  ; READ  IN  FIRST  16  ROWS 
CALL  CHECK! IER) 

TQG0LE«0  »  FLAG  TO  SHOW  BUFFER  WRAF -AROUND 

DO  9  1*1.7  i  COUNTER  FOR  4K  BUFFER  LOAOS 

CALL  RDBLK! INFIL. 16*1.  I NMAT < 1 . 1 7-16*T0CCLE) .  16.  IER) 

CALL  CHECK! IER)  »  READ  IN  NEXT  16  ROUS 

DO  3  KM.  16  i  COUNTER  FOR  OUTPUT  BUFFER  ROUS 

DO  2  JM,  8  >  ZERO  FILL  FIRST/LAST  8  COLUMNS 
OUTMAT  (J.K)-O 
OUTMAT!  JM 20,  K)-0 

2  CONTINUE 

DO  3  J-9, 120  t  ENHANCE  THIS  ROW 

CALL  E00ER1 ! J. K. TOGGLE. INMAT.  OUTMAT) 

3  CONTINUE 

CALL  WRBLK!OUTFIL.  16*1-8.  OUTMAT, 16.  IER)  i  WRITE  16  ROWS 
CALL  CHECK (IER)  i  OF  RESULTS 

IF (TOGGLE.  EO.  1)00  TO  4 

TOGGLE*!  i  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAO 
00  TO  3 

4  TOOOLE-O 
9  CONTINUE 

TYPE-EDGED  IMAGE  COMPLETE” 


RETURN 


SUBROUTINE  EDGER 1 < J, K. TOGGLE/  INMAT, OUTMAT) 

INTEGER  J, K. TOGGLE, Jl, J2, J3, Kt, K2, K3 

REAL  INMAT ( 12B,  32),  OUTMAT ( 128, 16) -  A,  B, C,  D,  E,  F,  C,  H 

IF (TOGGLE.  EG.  l)CO  TO  1 
K1-K4-7 
K2-K+8 
K3-K4-9 
GO  TO  2 

1  Kl«M0D<K+22. 32)4-1 
K2-MOD ( K+23.  32)4-1 
K3-MOO(K4-24.  32)4-1 

2  Jl-J-1 

J2-J 

J3-J4-1 

A«9*(  INMAT (J1.K1)4-1NMAT<  J2,  K1  >  4-INHAT  <J3.  Kl>>- 
4-  3*<  INMAT <J3,  K2 )  +  INMAT ( J3.  K3> 4-INMAT ( J2.  K3>4- 

4-  INMAT<J1.  K3)4-INMAT(J1,  K2)  ) 

B<“9*(  INMAT ( J2,  K1  )♦  INMAT (J3,  K1 )  +  INMAT ( J3,  K2>  )- 
4-  3*<  INMAT <  J3.  K3) + INMAT ( J2.  K3)*INMAT <  Jl,  K3>4- 

4-  INMAT(  Jl,  K2 )  4-INMAT  (Jl,  K1 )  ) 

C-9*( INMAT (J3,  K1  )4-INMAT(J3.  K2)4-INMAT<  J3,  K3>  )- 
4-  3*(  INMAT ( J2,  K3 )  +  INMAT <  Jl ,  K3) ♦INMAT ( Jl ,  K2)  4- 

4-  INMAT ( Jl,  K1 )♦ INHAT <J2»  K1 )  ) 

D-9*  ( I NMAT  ( J3.  K2 )  4- 1 NMAT  ( J3,  K3)+I  NMAT  ( J2,  K3 )  >  - 
4-  3*<INMAT(  Jl.  K3)4-INMAT(J1,  K2) 4-INMAT <J1,  Kl)4- 
4-  INMAT  ( J2,  K1 )  4- 1  NMAT  ( J3«  K1 )  ) 

E»9*(  INMAT ( J3,  K3)4-INMAT(  J2,  K3)  4-INMAT <  Jl ,  K3)  >- 
4-  3*(  INMAT(  Jl,  K2)4-INMAT(  Jl,  K1 ) 4-INMAT (J2,  Kl)4- 

4-  INMAT  ( J3.  K1 )  4- 1  NMAT  ( J3.  K2)  ) 

F»9«( INMAT (J2,  K3)4-INHAT(  Jl.  K3)4-INMAT(  Jl.  K2>  )- 
4-  3*(  I  NMAT  ( Jl.  XI  )t-INMAT(  J2.  K1 )  4-INMAT  <J3,  K1  )4- 

4-  INMAT(J3.  K2)4-INMAT(J3.  K3)> 

Q-9*(  INMAT(  Jl.  H3)4-INMAT(  Jl,  K2)4-INMAT(  Jl,  K1 )  >- 
♦  3*(  INMAT<  J2.  K1  )4-lNMAT(  J3.  K 1 )  ♦  I  NMAT  ( J3,  K214- 

4-  INMAT(J3.K3)4-INMAT(J2,K3)  > 

H-9*(  INMAT(  Jl,  K2)4-INMAT(  Jl.  K1  )4-INMAT(  J2,  K1  >  )- 
4-  3*(INMAT(J3.  K1)4-INMAT(J3.  K2)4-INMAT(J3.  K3)  + 

4-  INMAT(J2,  K3)4-INMAT(J1.K3)) 

OUTMAT ( J.  K)-AMAX1 (A.  B. C. D.  E,  F.  0.  H) 


RETURN 

END 


SUBROUTINE  C0NJUN(F1LE1.  F1LE2.  FILE3, ENTHR. EDTHR,  INMAT. OUTMAT) 
INTECER  FILE!. FILE2. FILE3 

REAL  ENTHR. EDTHR. INMAT < 128. 38) . OUTMAT ( 128, 16) 

DO  1  K— 1.  16  »  ZERO  FIRST/LAST  16  ROWS 

DO  1  J-l. 128 

OUTMAT !J. K)-0.  0 

1  CONTINUE 

CALL  UR9LK(FILE3.  0.  OUTMAT. 16. IER ) 

CALL  CHECK! IER) 

CALL  WRBLK( FILES.  112,  OUTMAT. 16. IER) 

CALL  CHECK! IER) 

DO  9  1-1.6  i  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  RDBLK!FILE1.  16*1,  INMAT! 1. 1).  16.  IER)  i  READ  16  ROWS 
CALL  CHECK! IER) 

CALL  R0BLK!FILE2. 16*1. INMAT! 1. 17), 16. IER)  ;  READ  16  ROWS 
CALL  CHECK! IER) 

DO  4  K-l. 16  i  COUNTER  FOR  OUTPUT  BUFFER  ROWS 

DO  3  J-17, 112  i  APPLY  CONJUNCTIVE  THRESHOLD  TEST 
IF! INMAT! J.  K).  QT.  ENTHR.  AND. 

♦  INMAT! J. K*16).  ©T.  EDTHR >CO  TO  2 

OUTMAT! J.  K >— O.  0 
00  TO  3 

2  OUTMAT! J.  K)  — 15.  0 

3  CONTINUE 

DO  4  J-l. 16  )  ZERO  THE  PERIMETER 

OUTMAT! J.K)-0 
OUTMAT !J+ 112. K)-0 

4  CONTINUE 

CALL  WRBLK!FILE3.  16*1.  OUTMAT,  16. IER)  i  WRITE  16  ROWS 
CALL  CHECK! IER) 

9  CONTINUE 

DO  6  K-l.  2  i  ZERO  THE  PERIMETER 

DO  6  J-l. 128 

OUTMAT! J.  l>-0.  O 
6  CONTINUE 

CALL  WRBLK!FIL£3.  16.  OUTMAT, 2,  IER) 

CALL  CHECK! XER) 

CALL  WRBUMFILE3.  110.  OUTMAT.  2.  IER) 

CALL  CHECK! IER) 

TYPE-CONJUNCTIVE  THRESHOLDING  DONE- 


RETURN 

END 


SUBROUTINE  CONNEC! INFIL.  OUTFIL.  INMAT. OUTMAT) 

INTEGER  INFIL, OUTFIL.  TOGGLE.  ICNT, JCNT 
REAL  INMAT ( 128. 32).  OUTMAT ( 128.  16) 

JCNT-1  j  COUNTER  FOR  TEST  LOOP  COUNTING 

1  ICNT-0  i  COUNTER  FOR  TEST  STABILITY 

CALL  ROBLK< INFIL. 0.  INMAT. 16.  IER)  ; READ  IN  FIRST  16  ROUS 
CALL  CHECK (IER) 

IF (JCNT. QT.  11INFIL-0UTFIL  > OPERATE  ON  SAME  FILE  AFTER  FIRST  LOOP 
TOGCLE-O  »  FLAG  TO  SHOW  BUFFER  WRAP-AROUND 

DO  9  1*1.7  i  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  ROBLK( INFIL.  16*1. INMAT ( 1. 17-16*T00CL£>. 16. IER) 

CALL  CHECK! IER)  ;  READ  NEXT  16  ROUS 

DO  3  K-l.  16  i  COUNTER  FOR  OUTPUT  BUFFER  ROUS 

DO  2  J-l.B  i  ZERO  FILL  FIRST/LAST  8  COLUMNS 
OUTMAT  (J.K)-O.  O 
OUTMAT! J+ 120.  K)*0.  0 

2  CONTINUE 

DO  3  J-9. 120  t  TEST  THIS  ROU  FOR  CONNECTIVITY 
CALL  C0NEC1 ( J. K, TOGGLE,  ICNT.  INMAT.  OUTMAT) 

3  CONTINUE 

CALL  URBLK( OUTFIL.  16*1-8.  OUTMAT,  16.  IER)  i  WRITE  16  ROUS 
CALL  CHECK! IER)  ;  OF  RESULTS 

IF(TOOOL£.  EO.  1)00  TO  4 

TOQQLE* 1  i  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 
GO  TO  9 

4  TOOOLE-O 
9  CONTINUE 

DO  6  K-l.  8  i  ZERO  FILL  FIRST/LAST  8  ROUS 
DO  6  J-l.  128 

OUTMAT  (J.K)-O.  0 
A  CONTINUE 

CALL  URBLK(0UTF1L, 0.  OUTMAT,  8,  IER) 

CALL  CHECK! IER) 

CALL  URBLK! OUTFIL, 120.  OUTMAT.  8,  IER) 

CALL  CHECK! IER) 

JCNT— JCNT*1 

IF! ICNT.  NE.  0)00  TO  1  »  CONNECTIVITY  TEST  HAS  NOT  STABILIZED 

TYPE"CONNECTIVITY  TEST  DONE  AT  LOOP  JCNT-1 

RETURN 

END 


0-9 


•  y 
■y 


SUBROUTINE  CONEC1 < J,  K,  TOOOLE.  ICNT.  INMAT, OUTMAT) 

INTEGER  J.  K. TOGGLE.  ICNT 

REAL  INMAT ( 128.  32) •  OUTMAT ( 128.  16) 

INTEOER  U1.U2.  U3.  K1.K2.K3.  NCNT 

IF ( TOOOLE.  EO.  1)00  TO  1 
Kl-K+7 
K2-K-8 
K3-K** 

00  TO  2 

1  Kl-M0D(K+22<  J32)+l 
K2-M0D<m-23.  32>fl 
K 3-MOD ( K+24,  32)4-1 

2  Jl-U-1 
J2— J 
J3—-J+1 

IF<INMAT<J2.  K2>.  EO.  0.  0)00  TO  4 
NCNT— 0 

IF( INNAT( J2<  K1 ) .  NE.  0.  0)NCNT-NCNT+1 
IF <  INMAT <J3,  K2).  NE.  0.  O ) NCNT-NCNT*- 1 
IF< INMAT <J2.  K3).  NE.  O.  0)NCNT-NCNT+1 
IF  < INMAT ( Jl.  K2) .  NE.  0.  O ) NCNT-NCNT+ 1 

IF (NCNT.  LT.  2)00  TO  3 
OUTMAT ( J.  K)— 19.  0 
00  TO  9 

3  OUTMAT ( J,  K ) -O 

ICNT— ICNT+1  ,  INCREMENT  COUNTER  TO  SHOW  CHANGE 
OO  TO  3 

4  OUTMAT < J.  H)-0.  0 

9  CONTINUE 

RETURN 

END 


PROGRAM  HAM21  (SPATIAL  CLUSTER  RECOGNITION  ALGORITHM) 

INTEGER  MAIN(7),  IFLNM(7)  .  F2(7) .  F3<  7)  i  MS(2)  •  SI (2) «  S2( 2) ,  S3(2> 
REAL  INMAT < 128.  32) .  QUTMAT ( 12S.  16) . ENPCT, EDPCT 

CALL  IOF( 1.  MAIN.  IFLNM.  F2.  F3,  MS. SI. S2. S3) 

CALL  OPENd.  IFLNM.  1.  IER 1 
CALL  CHECK (IER) 

CALL  OPEN <2.  "ENHC.  IR".2.  IER) 

CALL  CHECK (IER) 

CALL  OPEN (3.  "NTHR.  TH-.2.  IER) 

CALL  CHECK (IER) 

CALL  OPEN (4.  "NTHR.  IR“.  3.  IER) 

CALL  CHECK (IER) 

CALL  OP£N( 3< "EDGE.  IR-.2.  IER) 

CALL  CHECK (IER) 

CALL  OPEN (6. "DTHR.  TH". 2. IER) 

CALL  CHECK! IER) 

CALL  OPEN (7.  "DTHR.  IR".3. IER) 

CALL  CHECK! IER > 

CALL  OPEN (8.  “CTHR.  IR".2. IER) 

CALL  CHECK (IER) 

CALL  OPEN (9.  “CONN.  IR",2,  IER) 

CALL  CHECK! IER) 

ENHANCED  IMAGE  IS  FORMED  FROM  ORIGINAL  IMAGE 
CALL  NHANCE! 1.  2.  INMAT. OUTHAT) 

ENPCT-1.  0 
MASKSZ-17 

CALL  LTHRSH ( 2.  3.  4.  MASKSZ.  ENPCT.  INMAT,  OUTMAT) 

EDQED  IMAGE  IS  FORMED  FROM  ORIGINAL  IMAGE 
CALL  EOIMAO( 1,  3,  INMAT. OUTMAT) 

EDPCT- 1  0 

CALL  LTHRSH! 3.  6.  7,  MASKSZ.  EDPCT.  INMAT.  OUTMAT) 

EDOED  IMAGE  IS  REDUCED  BY  CONJUNCTIVE  THRESHOLDING 
CALL  LOCONJ ( 2.  3.  3.  6.  8.  INMAT.  OUTMAT) 

EDGED  IMAGE  IS  REDUCED  BY  CONNECTEDNESS  TEST 
CALL  CONNEC'B.  9.  INMAT,  OUTMAT) 

CALL  RESET 

8T0P"<7X7X7X7>HAH21  ‘ 

END 


SUBROUTINE  LTHRSH< INFILE.  OFILEi. 0FILE2. MASKS 2.  PCT.  INMAT,  OUTMAT) 

INTEGER  INFILE, OFILEI, 0FILE2, MASKSZ, TOGGLE, OFF 

REAL  PCT. INHAT! 128.  32) . OUTMAT < 128, 16) . MSUM,  NR SUM,  CHANGE 

DO  1  K-l,  8  i  FILL  FIRST  AND  LAST  8  ROWS  OF 

DO  1  J-l.  128  i  THRESHOLD  IMAGE  WITH  CONSTANT 

OUTMAT !  J.  K)«l.  0 

1  CONTINUE 

CALL  WRBLH < OFILEI, 0. OUTMAT, 8. IER) 

CALL  CHECK! IER) 

CALL  WRBLK(0FILE1, 120,  OUTMAT.  8.  IER) 

CALL  CHECK! IER) 

CALL  RDBLK ! INFILE, 0.  INMAT, 32,  IER)  {READ  IN  FIRST  32  ROWS 
CALL  CHECK! IER) 

OFF*! MASKSZ- 1 ) /2 
MSUM-O.  0 

DO  2  K— O. MASKSZ— 1  i  INITIALIZE  AVERAGING  VARIABLES 
DO  2  U-0, MASKSZ-1 

MSUM-MSUM+ 1 NMAT ! 9-OFF-t-J,  9-OFF ♦K ) 

2  CONTINUE 
CHANGE-0 

DO  3  J-0. MASKSZ-1 

CHANOE-CHANCE+INMAT ! 9— QFF+J.  9+OFF) 

3  CONTINUE 

NRSUM— M3UM— CHANCE 

TOOGLE-O  i  FLAO  TO  SHOW  BUFFER  WRAP-AROUND 

DO  4  K-l,  16  i  COMPUTE  THRESHOLD  FOR  FIRST  16  ROWS  !4K) 

CALL  THR0W1 !K.  MASKSZ. TOGOLE. NRSUM, PCT. INMAT.  OUTMAT) 

4  CONTINUE 

CALL  WR8LM0FILE1.8.  OUTMAT.  16.  IER)  »  WRITE  RESULTS 
CALL  CHECK! IER) 

TOGOLE- 1 

DO  7  1-2.7  i  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  RDBLK! INFILE. 16*1. INMAT! 1, 17-16*T00GLE> .  16,  IER) 

CALL  CHECK! IER)  i  READ  IN  NEXT  16  ROUS 

DO  9  K-l. 16  >  COUNTER  FOR  OUTPUT  BUFFER  R0W8 

CALL  THROW 1 !K.  MASKSZ. TOGGLE,  NRSUM.  PCT.  INMAT.  OUTMAT) 

3  CONTINUE 

CALL  WRBLK! OF I LEI. 16*1-8.  OUTMAT. 16, IER)  i  WRITE  16  ROUS 
CALL  CHECK  HER) 

IF ! TOGGLE.  EQ.  1)60  TO  6 

TOGOLE- 1  i  TOGOLE  THE  BUFFER  WRAP-AROUND  FLAO 
GO  TO  7 

6  TOOGLE-O 

7  CONTINUE 

DO  9  1-0,7  i  FORM  THRESHOLDED  IMAOE 

CALL  RDBUM0F1LE1,  16*1.  INMAT,  16.  IER) 

CALL  CHECK! IER) 

CALL  RDBLK! INFILE,  16*1.  OUTMAT, 16. IER) 

CALL  CHECK! IER) 


PACE 

00  8  J-l. 128 

IF! OUTMAT <  J.  K ) .  LE.  INHAT! J. K> )OUTMAT!J. K)»0.  0 

8  CONTINUE 

CALL  WRBLK ( OF I LE2<  16*1. OUTMAT. 16. IER ) 

CALL  CHECK! IER) 

9  CONTINUE 

TYPE " THRESHOLD  MATRIX  COMPLETE" 

RETURN 

END 


SUBROUTINE  THROW1 <K, MASKSZ. TOGGLE. NR SUM. PCT. INMAT. OUTMAT) 


INTEGER  K, MASKSZ. TOGGLE,  OFF 

REAL  NRSUM. PCT,  INMAT< 129,  32 ) ,  OUTMAT < 128.  16) 
REAL  FCOLSM.  CHANGE.  ASUM,  AUG.  DSUM,  DEV 


IF (TOGGLE.  EG.  l)GO  TO  1 
Kl»7 
00  TO  2 

1  Kl-23 

2  CONTINUE 

OFF- ( MASKSZ- l)/2 


INITIALIZE  INDEXING  VARIABLES 


DO  3  J-1.8  i  FILL  FIRST/LAST  EIGHT  PIXELS  WITH  CONSTANT 
OUTMAT  <  J,  K )  — 1.  0 
0UTMAT(J*120.  K)-l.  0 

3  CONTINUE 

CHANGE-0.  0  i  FIND  NEW  MATRIX  SUM 

DO  4  1-0. MASKSZ-1 

CHANCE -CHANGE+INMAT <9-0FF*I»  M0D(K*0FF+K1.  32>*1 ) 

4  CONTINUE 
ASUM-NRSUM+C MANGE 

CHANOE-O.  0  i  FIND  NEW  NEXT  ROW  SUM 

DO  3  1-0. MASKSZ-1 

CHANGE-CHANGE* INMAT ( 9-OFF* I, M0D(K-0FF*K1,  32)*1 ) 

3  CONTINUE 

NRSUM— ASUM— CHANCE 

FCOLSM— O.  0  i  INITIALIZE  FCOLSM 

DO  6  1-0. MASKSZ-1 

FCOLSM— FCOLSM* INMAT <  9*0FF»  M0D(K-0FF*I*K1,  32) *1 ) 

6  CONTINUE 

DO  10  J— 9,  120  i  COLUMN  COUNT 

CHANGE-0.  O  ;  FIND  NEW  MATRIX  SUM  AND  AVERAGE 
DO  7  I -0, MASKSZ-1 

CHANGE-CHANGE* I NMAT ( J*OFF ,  M0D<K-0FF+I*K1.  32>*1> 

7  CONTINUE 

ASUM— ASUM+CHANCE -FCOLSM 
AVO— ASUM/ MASKSZ  *■*2 

FCOLSM— 0.  O  i  FIND  NEW  FIRST  COLUMN  SUM 

DO  8  1-0. MASKSZ-1 

FCOLSM— FCOLSM* I NMAT  <  J-OFF .  MOD ( K-OFF* I *K 1 .  32)*1 ) 

8  CONTINUE 

DSUM— O.  0  »  FIND  NEW  DEVIATION 

DO  9  I 1-0, MASKSZ-1 
DO  9  12-0, MASKSZ-1 

DSUM— DSUM* < I NMAT  <  J-OFF* 1 1 . MOD ( K -OFF* 1 2*K 1 .  32 ) * 1 ) -AVO ) **2 

9  CONTINUE 
DEV-SGRT(DSUM/MASKSZ**2) 

OUTMAT (J.  K)-AVG*PCT*DEV 

10  CONTINUE 

RETURN 

END 
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SUBROUTINE  LOCONJ(FILEl*  FILE2, FILES. FILE4, FILES,  INMAT, OUTMAT) 

INTEGER  FILE1, FILE2,  FILE3. FILE4. FILES 
REAL  INMAT < 12S, 32),  OUTMAT< 128,  16) 

00  1  K-l. 16  i  ZERO  FIRST/LAST  16  ROWS 
00  1  J-l. 128 

OUTMAT !J. K)-0.  0 
1  CONTINUE 

CALL  WRBLK (FILES, 0,  OUTMAT. 16,  I£R> 

CALL  CHECK (1ER) 

CALL  WR8LK!FILE9.  112.  OUTMAT,  16.  IER) 

CALL  CHECK! IER ) 

DO  S  11-1.6  i  COUNTER  FOR  4K  BUFFER  OUTPUT 

00  4  12-0. 1  s  COUNTER  FOR  2K  BUFFER  INPUT 

CALL  RDBLK!FILE1.  16*11+8*12.  INMAT(1,  1), 8.  IER)  ;  READ  8  ROWS 
CALL  CHECK! IER) 

CALL  RDBLK(F1LE2,  16*11+8*12,  INMAT! 1. 9) . 8.  IER)  i  READ  8  ROWS 
CALL  CHECK! IER) 

CALL  RDBLK!F!L£3>  16*11+8*12.  INMAT! 1.  17). 8. IER)  s  READ  8  ROWS 
CALL  CHECK! IER) 

CALL  RDBLK!F1LE4, 16*1 1+8*I2.  INMAT  1 1 , 29) . 8. IER)  >  READ  8  ROWS 
CALL  CHECK! IER) 

DO  4  K-l. 8  i  COUNTER  FOR  OUTPUT  BUFFER  ROWS 

DO  3  J-17.  112  i  APPLY  CONJUNCTIVE  THRESHOLD  TEST 

IF!  INMAT!  J,K).  OT.  INMAT!  J.  K«-B> .  AND. 

♦  INMAT! J, K*16>.  OT.  INMAT! J, K+24> >00  TO  2 

OUTMAT! J.  K*B»I2)-0.  0 
00  TO  3 

OUTMAT!  J.  K«-B*I2)-19.  0 
CONTINUE 

DO  4  >1.  16  i  ZERO  THE  PERIMETER 
OUTMAT! J.  K*8*I2)-0.  0 
OUTMAT! J*1 12.  M-8»I2)-0.  O 


4  CONTINUE 

CALL  URBLK!FILE9,  16*11.  OUTMAT.  16.  IER)  i  WRITE  16  ROWS 
CALL  CHECK! IER) 

9  CONTINUE 

DO  6  K-l. 2  i  ZERO  THE  PERIMETER 

DO  6  J-l. 128 

OUTMAT ! J, K ) «0.  0 
6  CONTINUE 

CALL  WRBLK!FILE9. 16. OUTMAT. 2.  IER) 

CALL  CHECK! IER) 

CALL  WRBLK! FILES. 1 10.  OUTMAT,  2,  IER) 

CALL  CHECK! IER) 

TYPE "LOCAL  CONJUNCTIVE  THRESHOLDING  DONE" 

RETURN 

END 
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PROGRAM  HAM31  (SPATIAL  CLUSTER  RECOGNITION  ALGORITHM) 


INTEGER  MAIN<7>. IFLNM<7>. F2(7>. F3!7>. MS! 2). SI <2). S2(2).  S3!2> 
REAL  INMAT ( 128. 32)  •  OUTMAT ( 12B>  16 > . PCT,  DMIN, DMAX 

CALL  IOF( 1. MAIN. IFLNM. F2. F3,  MS.  SI .  S2,  S3 > 

CALL  OPEN ( 1 < IFLNM. 1. IER) 

CALL  CHECK! IER) 

CALL  OPEN ( 2.  "ENHC.  IR-.2.  IER) 

CALL  CHECK (IER) 

CALL  OPEN! 3.  "NTHR.  TH".2.  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 4. "NTHR.  IR",  3. IER) 

CALL  CHECK! IER) 

CALL  OPEN! 9.  "EDGE.  IR",2,  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 6.  “DTHR.  TH".  2.  IER) 

CALL  CHECK! IER) 

CALL  0PEN!7»  "DTHR.  IR". 3.  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 8. "CTHR.  IR".  2.  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 9. "CONN.  IR".  2.  IER) 

CALL  CHECK! IER) 

C  ENHANCED  IMAGE  IS  FORMED  FROM  ORIGINAL  IMAGE 

CALL  NHANCE! 1.  2.  INMAT.  OUTMAT) 

MASKSZ-7 
PCT-1.  0 
DMIN-IOO.  0 
DMAX- 190.  0 

CALL  L0THR8 ! 2. 3.  4.  MASKSZ. PCT.  DMIN. DMAX.  INMAT. OUTMAT) 

C  EDOED  IMAGE  IS  FORMED  FROM  ORIOINAL  IMAGE 


CALL  EOIMAO! 1.  9. INMAT.  OUTMAT) 

DMIN— O.  O 
DMAX— 10.  0 

CALL  L0THRS!9. 6. 7.  MASKSZ. PCT.  DMIN. DMAX. INMAT. OUTMAT) 
EDOED  IMAGE  IS  REDUCED  BY  CONJUNCTIVE  THRESHOLD I NO 
CALL  LQCONJ ! 2. 3.  9.  A.  8.  INMAT.  OUTMAT) 

EDGED  IMAGE  IS  REDUCED  BY  CONNECTEDNESS  TEST 


CALL  CONNEC ! 8. 9. INMAT.  OUTMAT) 
CALL  RESET 

STOP  "  <7><7X7><7>MAM3 1 " 

END 
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SUBROUTINE  LOTHRS< INFILE. OFILE1.  OFILE2.  MASKSZ.  PCT,  DMIN.  DMAX. 

♦  INMAT. OUTMAT) 

INTEGER  INFILE. OFILE1. 0FILE2.  MASKSZ. TOGOLE.  OFF 

REAL  PCT.  DMIN.  DMAX.  INMAT ! 128.  32 > .  QUTMAT( 128.  16).  MSUM.  NRSUM. CHANGE 

DO  1  K-l,  8  i  FILL  FIRST  AND  LAST  8  ROWS  OF 

DO  1  J-l. 128  i  THRESHOLD  IMAGE  WITH  CONSTANT 

OUTMAT! J.  K>-1.  0 
1  CONTINUE 

CALL  WRBLK<OFILEl.  0.  QUTMAT.  8.  IER) 

CALL  CHECK < IER) 

CALL  WRBLK< 0FILE1 .  120.  OUTMAT,  8.  IER) 

CALL  CHECK (IER) 

CALL  RDBLM INFILE. 0, INMAT. 32.  IER)  .READ  IN  FIRST  32  ROWS 
CALL  CHECK! IER) 


OFF* <  MASKSZ -1 ) /2 
MSUM-O.  0 

DO  2  K*0.  MASKSZ-1  »  INITIALIZE  AVERAGING  VARIABLES 
DO  2  J-0, MASKSZ-1 

MSUH-MSUN+ 1 NMAT <  9-QFF+ J.  9-OFF+K > 

2  CONTINUE 
CHANGE-0 

DO  3  J-0.  MASKSZ-1 

CHANOE-CHANCE+INMAT  <  9-OFF+J,  9-t-OFF  > 

3  CONTINUE 

NR SUM-MSUM— CHANGE 

TOGOLE— O  i  FLAG  TO  SHOW  BUFFER  WRAP-AROUND 


DO  4  K-l,  16  i  COMPUTE  THRESHOLD  FOR  FIRST  16  ROWS  <4K) 

CALL  THR0W2(K. MASKSZ.  TOGOLE. NRSUM. PCT.  DMIN.  DMAX.  INMAT.  OUTMAT) 
4  CONTINUE 

CALL  WRBLK(0FILE1. 8.  OUTMAT, 16.  IER)  ;  WRITE  RESULTS 
CALL  CHECK (IER) 

TOGGLE- 1 

DO  7  1-2. 7  i  COUNTER  FOR  4K  BUFFER  LOADS 

CALL  RDBLK< INFILE.  16*1. INMAT < 1. 17-16*T000LE>,  16.  IER) 

CALL  CHECK! IER)  >  READ  IN  NEXT  16  ROWS 


DO  9  K-l. 16  i  COUNTER  FOR  OUTPUT  BUFFER  ROWS 

CALL  THROWS !K.  MASKSZ.  TOGGLE.  NRSUM.  PCT,  DMIN.  DMAX. 

♦  INMAT, OUTMAT) 

9  CONTINUE 

CALL  WRBLK! OF I LEI.  16*1-8.  OUTMAT. 16.  IER)  I  WRITE  16  ROWS 
CALL  CHECK! IER) 

IF (TOGOLE.  EO.  1)00  TO  6 

TOGGLE-1  i  TOGGLE  THE  BUFFER  WRAP-AROUND  FLAG 

GO  TO  7 

6  TOGGLE -O 

7  CONTINUE 

DO  9  1-0.7  <  FORM  THRESHOLDED  IMAGE 

CALL  RDBLK(0FILE1.  16*1, INMAT, 16. IER) 

CALL  CHECK! IER) 

CALL  RDBLK! INFILE.  16*1. OUTMAT, 16, IER) 
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CALL  CHECK  <  IER ) 


■y  v.% 
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DO  8  H-l. 16 
DO  8  J-l.  128 

IF<OUTMAT<J<  K).  LE.  INMAT< J. K> )OUTMAT< J.  K>-0  0 

8  CONTINUE 

CALL  WRBU<<0FILE2< 16*I.0UTMAT. 16. IER) 

CALL  CHECK < IER) 

9  CONTINUE 

TYPE" THRESHOLD  MATRIX  COMPLETE" 

RETURN 

END 
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SUBROUTINE  THR0W2CK,  MASKSZ,  TOGGLE,  NR SUM,  PCT,  DM1N.  DMAX, 

♦  INMAT, OUTMAT ) 

INTEGER  K.  MASKS Z,  TOGGLE.  OFF 

REAL  NRSUM, PCT,  DMIN.  DMAX.  INMAT< 128.  32),  OUTMAT< 128.  16) 

REAL  FCOLSM, CHANGE,  A SUM.  AVO.  DSUM,  DEV 

IF ( TOGGLE.  EG.  l)GO  TO  1  i  INITIALIZE  INDEXING  VARIABLES 
Kl-7 
00  TO  2 

1  Kl-23 

2  CONTINUE 

OFF*  <  MASKSZ-l )  /2 

DO  3  J-1.8  i  FILL  FIRST/LAST  EIOHT  PIXELS  WITH  CONSTANT 
OUTMAT <J,  K>-1. 0 
OUTMAT ( J+120,  K)— 1.  0 

3  CONTINUE 

CHANGE-0.0  i  FIND  NEW  MATRIX  SUM 

DO  4  1-0, MASKSZ-l 

CHANGE -CHANGE+INMAT (9-OFF+I,  M0D(K-H3FF+K1.  32H-1 ) 

4  CONTINUE 

ASUM— NRSUH+CHANOE 

CHANCE-0  0  i  FIND  NEW  NEXT  ROW  SUN 

DO  9  1-0.  MASKSZ-l 

CHANOE-CHANGE+INMAT(9-OFF+I,  MOD < K-OFF+K1 , 32)  >1) 

9  CONTINUE 

NRSUH— ASUH— CHANOE 

FCOLSM— 0.  O  s  INITIALIZE  FCOLSM 

DO  6  1-0.  MASKSZ-l 

FCOLSM— FCQLSH+  INMAT  < 9 ♦OFF,  MOD  <  K-OFF+ 1  ♦*  1 . 32  )  ♦  1 ) 

6  CONTINUE 

DO  10  J— 9,  120  i  COLUMN  COUNT 

CHANOE— O.  0  )  FIND  NEW  MATRIX  SUM  AND  AVERAGE 

DO  7  1-0. MASKSZ-l 

CHANOE-CHANOE+ INMAT <  J+OFF,  MOD(K-OFF+I^KI.  321  +  1 ) 

7  CONTINUE 
ASUM-ASUM+CHANOE-FCOLSM 
AVO-ASUM/ MASKS Z**2 

FCOLSM— O.  0  I  FIND  NEW  FIRST  COLUMN  SUM 

DO  8  1-0, MASKSZ-l 

FCOLSM— FCOLSM^ INMAT <  J-OFF,  MOD< K-OFF+I+Kl .  32>  +  l ) 

8  CONTINUE 

DSUM-O.  0  i  FIND  NEW  DEVIATION 
DO  9  I 1-0, MASK8Z-1 
DO  9  12-0,  MASKSZ-l 

DSUM— DSUM+<INMAT<U-OFF^Il,  MOO  <  K-OFF>  1 2-*K  1 .  32 >  + 1) -AVG ) **2 

9  CONTINUE 

DEV— SORT (DSUM/MASKSZ **2) 

D6V-D£V#PCT  i  MODIFY  DEVIATION 

IF(DEV.  LT.  DMIN) DEV-DMIN 
IF (DEV.  OT. DMAX ) DEV- DMAX 
OUTMAT (J, K) -A VO ♦DEV 

10  CONTINUE 

RETURN 

END 
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PROGRAM  HAM61  (SPATIAL  CLUSTER  RECOGNITION  ALGORITHM) 

INTEGER  MAIN( 7 ) «  IFLNM< 7) » F2<  7) » F3< 7) . MS ( 2) • SI ( 2) « S2(2) >  S3 ( 2 ) 
REAL  INMAT ( 128.  32 )  •  ENHMAT (128.  16). PCT, DMIN.  DMAX 

CALL  IOFll.MAIN,  IFLNH.  F2.  F3.  MS.  S1.S2.  S3) 

CALL  OPEN< 1. IFLNM. 1.  IER) 

CALL  CHECK! IER) 

CALL  OPEN <2.  “ENHC.  IR",2.  IER) 

CALL  CHECK! IER ) 

CALL  OPEN <3.  "NTHR.  IR", 2.  IER) 

CALL  CHECK! IER) 

CALL  OPEN! 4.  -CONN.  IR-.2.  IER) 

CALL  CHECK! IER) 

ENHANCED  IMAGE  FORMED  FROM  ORIGINAL  IMAGE 

CALL  NHANCE! 1. 2. INMAT.  ENHMAT) 

LOCAL  THRESHOLDS  COMPUTED  FOR  ENHANCED  IMAGE 
AND  ENHANCED  IMAGE  LOCAL  THRESHOLDEO 

MASKSZ-7 
PCT-O.  79 
DMIN-IOO.  O 
DMAX *1000000.  0 

CALL  NTHRSH ! 2.  3. MASKSZ. PCT. DMIN. DMAX.  INMAT. ENHMAT) 
THRESHOLDEO  IMAGE  IS  REDUCED  BY  CONNECTEDNESS  TEST 
CALL  C0NNEC!3< ♦. INMAT.  ENHMAT) 

CALL  RESET 

STOP “ <7><7><7><7>HAM6 1 " 

END 


3U3SGUTINE  NTi-.a3m£NhFIi_.  NIRFIl,  MASKSZ,  PCT,  DMIN,  DMAX,  InMAT.  ENHMAT 


INTEGER  ENHFIL.  NIRFIL,  MASKSZ 

REAL  PCT.  DMIN.  DMAX,  INMAT < 128. 32) . ENHMAT ( 128.  16) 

INTEGER  OFF, TOGGLE 
REAL  MSUM, NR  MSUM,  CHANGE 

CALL  RDBLK (ENHFIL,  8.  INMAT,  32,  IER)  ;  READ  IN  FIRST  32  ROUS 
CALL  CHECK (IER)  ;  OF  ENHANCED  IMAGE 

0FF-(MASKSZ-l)/2 
MSUM-O.  0 

DO  1  K-O,  MASKSZ-1  i  INITIALIZE  AVERAGING  VARIABLES 
DO  1  J-O. MASKSZ-1 

MSUM-MSUM+ 1 NMAT  < 1 7-OFF +J. 9-OFF-HO 

1  CONTINUE 
CHANGE *0.  0 

DO  2  J-O. MASKSZ-1 

CHANOE-CHANCE*INMAT( 17-OFFfJ,  9+OFF) 

2  CONTINUE 
NRMSUM— MSUM-CHANCE 

TOCOLE-O  )  FLAG  TO  SHOW  BUFFER  URAP-AROUND 

DO  3  K-l,  16  i  THRESHOLD  FIRST  16  ROUS 

CALL  NTR0U1 (K,  MASKSZ.  TOGGLE,  NR MSUM,  PCT. DMIN. DMAX.  INMAT,  ENHMAT ) 

3  CONTINUE 

CALL  URBLK (NIRFIL.  16.  ENHMAT. 16.  IER)  )  UR I TE  RESULTS 
CALL  CHECK (IER) 

TOGGLE- 1 

DO  6  1-2. 6  >  THRESHOLD  ENHANCED  IMAGE 

CALL  RDBLK (ENHFIL.  l6«l+8.  INMAT ( 1.  17-16*TQQGLE> . 16.  IER) 

CALL  CHECK ( IER )  <  READ  IN  NEXT  16  ROUS 

DO  4  K-l. 16  i  COUNTER  FOR  OUTPUT  BUFFER  ROUS 

CALL  NTR0H1 (K.  MASKSZ.  TOGGLE.  NR MSUM. PCT,  DMIN. DMAX, 

♦  INMAT.  ENHMAT) 

4  CONTINUE 

CALL  URBLK(NIRFIL.  16*1.  ENHMAT,  16.  IER)  i  WRITE  16  ROUS 
CALL  CHECK (IER) 

IF( TOGGLE.  Ed.  1)00  TO  9 

TOGGLE- 1  i  TOGGLE  THE  BUFFER  URAP-AROUND  FLAG 
GO  TO  6 

9  TOGOLE-O 

6  CONTINUE 

DO  7  K-l.  16  i  ZERO  FILL  FIRST/LAST  18  ROWS 
DO  7  J-l. 128 

ENHMAT (J.K)-O.  0 

7  CONTINUE 

CALL  URBLK ( NIRFIL.  0,  ENHMAT, 16.  IER) 

CALL  CHECK! IER) 

CALL  URBLK (NIRFIL.  16.  ENHMAT.  2,  IER  > 

CALL  CHECK! IER) 

CALL  URBLK(NIRFIL.  1 10.  ENHMAT, 2, IER) 

CALL  CHECK! IER) 

CALL  URBLK (NIRFIL,  112.  ENHMAT. 16, IER) 

CALL  CHECK! IER) 

TYPE " THRESHOLDED  ENHANCED  IMAOE  COMPLETE" 

RETURN 

END 
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SUBROUTINE  NTROWl <K, MASKSZ.  TOGGLE- NRSUM, PCT,  DMIN.  DMAX, 

*  INMAT, OUTMAT) 

INTEGER  K. MASKSZ, TOGGLE,  OFF 

REAL  NRSUM,  PCT, DMIN.  DMAX.  INMAT ( 128, 32), OUTHAT ( 128.  16) 
REAL  FCOLSM, CHANGE, ASUM, AUG,  DSUM.  DEU 

IF (TOGGLE.  EG.  1)60  TO  1  j  INITIALIZE  INDEXING  VARIABLES 
Kl-7 
GO  TO  2 

1  Kl-23 

2  CONTINUE 
0FF-(MASKSZ-l>/2 

DO  3  J-l, 16  i  FILL  FIRST/LAST  16  PIXELS  WITH  ZERO 
OUTMAT ( J. K>-0.  0 
OUTMAT  <  J* 112,  K>-0.  0 

3  CONTINUE 

CHANOE-O.  0  i  FIND  NEW  MATRIX  SUM 

DO  4  1-0, MASKSZ-1 

CHANCE-CHANQE+INMAT ( 17-0FF*I,  M0D(K*0FF*K1,  32>+l> 

4  CONTINUE 
ASUM-NRSUM+CHANOE 

CHANGE-0.  0  j  FIND  NEW  NEXT  ROW  SUM 

DO  9  1*0.  MASKSZ-1 

CHANGE-CHANGE-* INMAT <  17-OFF* I,  MOD < K-OFF+K 1 ,  32)*1 ) 

9  CONTINUE 

NRSUM— ASUM— CHANGE 

FCOLSM— O.  O  )  INITIALIZE  FCOLSM 

DO  6  1-0. MASKSZ-1 

FCOLSM— FCOLSH*INMAT <  17+OFF,  M0D<K-0FF*I*K1,  32)*1 ) 

6  CONTINUE 


DO  10  J— 17, 112 


i  COLUMN  COUNT 


CHANGE-0.0  J  FIND  NEW  MATRIX  SUM  AND  AVERAGE 
DO  7  1-0. MASKSZ-1 

CHANGE -CHANGE* INMAT ( J+OFF, M0D(K-0FF*I+K1,  32>*1) 

CONTINUE 

ASUM— ASUM+CHANGE-FCOLSM 
AVO— ASUM/MASKSZ**2 

FCOLSM— O.  O  )  FIND  NEW  FIRST  COLUMN  SUM 

DO  8  1-0.  MASKSZ-1 

FCOLSM— FCOLSM* I NMAT <  J-OFF, MOD  <  K-OFF* I *K 1 , 32 ) ♦ 1 ) 

CONTINUE 

DSUM— 0.  0  i  FIND  NEW  DEVIATION 

DO  9  I 1-0, MASKSZ-1 
DO  9  12-0. MASKSZ-1 

DSUM-DSUM* I ABS ( I NMAT ( J-OFF* 1 1 .  MOD ( K -OFF* 1 2*K 1 .  32 ) * 1 ) -AVG  ) 
CONTINUE 

DEV— D6UM/MASKSZ**2 

DEV-D€V*PCT  i  MODIFY  DEVIATION 

IF(DEV.  LT.  DMIN) DEV— DMIN 
IF (DEV.  OT.  DMAX > DEV— DMAX 

OUTMAT (J.  K ) -O.  0  ,  THRESHOLD  RULE 

IF( I NMAT (J, MOD ( K*K 1 . 32>*1 )  QT.  (AVG*DEV> ) OUTMAT ( J,  K ) - 1 9.  0 


10  CONTINUE 
RETURN 


ill 
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ABSTRACT 

Edge  extraction  Is  an  Image  processing  technique  for  defining 
the  edge  information  in  an  Image.  This  effort  researches  different 
edging  processes  as  applied  to  preprocessing  for  two  pattern 
recognition  processes.  The  first  one  Is  a  cross-correlation  method 
to  find  a  target  given  that  the  target  has  a  known  size,  orientation, 
and  aspect.  Correlation  Is  performed  In  the  spatial  frequency  domain 
with  two-dimensional  fast  Fourier  transforms  of  the  searched  edge 
Image  and  a  hand  drawn  edge  template  to  correct  for  translation  only. 

The  second  pattern  recognition  process  researched  also  uses 
edging  as  one  step  of  a  purely  spatial  domain  algorithm.  The 
approach  locates  targets  in  Infrared  Images  that  can  be  described 
a  "hot"  clusters.  A  cluster  recognition  algorithm  by  Hamadanl  Is 
Implemented  and  altered  for  testing  of  local  thresholding  and 
thresholding  rules.  The  algorithm  Is  shown  to  be  effective  on 
real  Infrared  Images  provided  by  the  thesis  sponsor,  the  Air  Force 
Armament  Laboratory  at  Eglln  AFB,  Florida. 


